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Emile BOREL 
(7 janvier 1871 - 3 février 1956) 


En la personne d'Emile Borel, l'Institut de Statistique de l'Université 
de Paris a perdu le 3 février dernier le président de son Conseil 
d'Administration et les Publications de l'Institut, un des membres les plus 
éminents de son Comité Directeur. En même temps disparaissait une 
personnalité scientifique de premier plan dont l'influence à marqué 
profondément la production mathématique française du début de ce siècle. 
Emile Borel était l'homme de la l'réalité mathématique!'', Toute l'oeuvre 
du savant porte le signe de l'intérêt qu'il portait au caractère constructif 
d'une solution. Une fonction, un ensemble, c'était, pour lui, quelque chose 
que i'on doit pouvoir déterminer en un nombre fini de mots ou d'opérations, 


M. Arnaud Denjoy a écrit dans un commentaire publié dans le volume 
jubilaire de l'illustre mathématicien : ''Je doute qu'un seul analyste ait 
autant que M. Borel, ni même à un degré comparable ausien, le sens de la 
réalité numérique exprimée par les symboles mathématiques!'!, 


On connait la discussion célèbre connue sous le nom de querelle de 
l'axiome du choix à laquelle participèrent Baire, Lebesgue et M. Hadamard, 
Dans la circonstance, M. Borel fut le défenseur de la position extrême : 
L'existence en mathématique doit être réservée à ce qu'un mathématicien 
peut construire et peut décrire complètement à un autre mathématicien. 
Il semble d'ailleurs que dans ce domaine sa pensée ait un peu évolué. Un 
de ses derniers ouvrages, ''Les Nombres inaccessibles'', parait l'indiquer. 
Il était naturel qu'un esprit ayant ces dispositions, tourne ses puissantes 
facultés d'analyse vers cette science de l'application qu'est la Statistique. 
Un mémoire universellement connu : ‘Les probabilités dénombrables et 
leurs applications arithmétiques'' publié aux Rendiconti de Circolo 
Matematico de Palerme, avait marqué son passage de la théorie des 
fonctions oùils'étaitillustré en généralisantetapprofondissant les résultats 
d'Emile Picard et en conquérant la notion de fonction monogène, au Calcul 
des Probabilités.Ce mémoire attache ra le nom d'Emile Borel à la définition 
de la convergence presque certaine, conséquence toute naturelle de ses 
travaux sur la mesure des ensembles, généralisation que le grand analyste 
a voulu aussi simple et naturelle que possible de mesure des longueurs, 
Ces recherches l'avaient conduit à dégager deux sortes d'ensembles : 
les ensembles mesurables B (au sens de Borel) qui constituent ce que l'on 
appelle maintenant en calcul des probabilités le corps de Borel, et les 
ensembles de mesure nulle, La réunion de ces deux classes constitue les 
ensembles mesurables au sens de Lebesgue : seuls peuventêtre utilisables, 
d'une manière générale, en calcul des probabilités, les ensembles mesu- 
rables B.Il est saisissant de remarquer combien, dans une étude aussi vaste 
et aussi profonde, Emile Borel a lancé d'idées neuves sans surcharger le 
langage mathématique de mots nouveaux, 


Cette convergence presque certaine, distincte de la convergence en 
probabilité connue depuis Bernoulli, M. Fréchet l'a qualifiée dans son 
traité, ‘'d'acquisition la plus importante dans ce domaine depuis les travaux 
de Laplace'', Elle a entrainé tous les travaux de Kolmogoroff, P. Lévy, 
Robbins sur les différentes formes de loi des grands nombres que l'on peut 
en déduire. 


Au même moment, Emile Borel s'intéressait à la statistique appliquée 
et en 1909, à la demande d'Arthur Fontaine, alors directeur au Ministère 
du Travail, il participait à la session de Paris de l'Institut International de 
Statistique en présentant un travail sur l'emploi de la Méthode différen- 
tielle pour lacomparaisondes statistiques .Emile Borel devait rester fidèle 
toute sa vie aux réunions de ce grand organisme international : il a participé 
aux dernières réunions de New-Delhi, Rome, Rio de Janeiro. C'est la 29€ 
Session de l'Institut International dans cette dernière ville quia recueilli 
sa dernière communication sur l'étude statistique du vocabulaire des 
écrivains. L'Institut International devait en 1949 lui rendre hommage enle 
faisant figurer parmi ses membres honoraires : douze membres seulement 
ont à l'heure actuelle cette qualité et parmi eux les anciens présidents de 
République : H. Hoover des Etats-Unis et Einaudi d'Italie, ainsi que 
Sir Ronald Fisher, 


L'enseignement de la statistique en France avait été l'objet de ses 
préoccupations. En 1922, il fut un des fondateurs et un des premiers 
professeurs de l'Institut de Statistique de l'Université de Paris, l'année 
même où la Société de Statistique de Paris l'élisait comme Président, 
Il devait y conserver son enseignement jusqu'à son élection en 1924 
comme député de l'Aveyron. Il n'interrompit son cours de Calcul des 
Probabilités à la Faculté des Sciences de Paris que pendant l'année où 
son collègue, confrère et ami Paul Painlevé lui confia le Ministère de la 
Marine en 1925, C'est d'ailleurs pendant l'exercice de cette haute charge 
qu'il écrivit le fascicule de son traité relatif aux applications du Calcul 
des Probabilités, à l'Arithmétique, avec l'aide de M. Paul Dubreil 
maintenant professeur d'Arithmétique et de Théorie des Nombres à la 
Sorbonne, L'ensemble de tous les volumes de ce traité restera une 
lsomme'' dont la forme est aussi parfaite que dense est le fond, IL est 
certain que longtemps encore les jeunes probabilistes viendront s'inspirer 
de ces leçons dont M. Jean Ville, qui a rédigé celles qui sont relatives à 
la théorie des jeux, a dit dans son discours jubilaire de janvier 1940 


Ce qui nous attira dès l'abord vers les ouvrages de M. Borel, je le 
dirai en toute simplicité : ils se lisaient sans effort, et ils nous portaient, 
d'un élan rapide de notre savoir scolaire jusqu'aux travaux originaux et 
aux traités spécialisés!!, 


Un esprit aussi vaste ne saurait se résumer en un nombre fini de 
mots. En concluant cependant, il faut rappeler le profond patriotisme 
d'Emile Borel. Agé en 1915 de 44 ans, il avait demandé à commander 
une batterie d'artillerie, Sa conduite au Chemin des Dames,où il avait 
replié son unité après avoir épuisé toutes ses munitions et sans perdre 
un homme lui avait valu la Croix de guerre, En 1940,malgréles rigueurs 
de l'occupation qui allaient entrainer pour lui à plus de 70 ans l'arrestation 
et la détention, il avait regagné son poste à Paris pour reprendre son 
cours, 


La rédaction des publications de l'Institut de Statistique de Paris 
s'incline très respectueusement devant le souvenir de cette grande figure 
et présente à Madame Emile Borel l'hommage de ses condoléances. 


DISTANCE - DISCRIMINATION 
ET RÉSUMÉ EXHAUSTIF 


par 


B. P. ADHIKARI et D. D. JOSHI 


INTRODUCTION 


A. RÉSUMÉ HISTORIQUE 


C'est avec Karl Pearson (voir [20]) qu'est commencée l'étude d'une 
distance entre deux univers statistiques. Pearson a proposé une distance 
entre deux ensembles de mesures anthropométriques, et il l'a appelée le 
l'coefficient de similarité raciale''. Soit 


(Cxi;)) i=l, ..., p; j=l, ...,n: 
et 
((yx)) MES DURE INE A Unr 


deux tels ensembles de mesures où les mesures du caractère i viennent de 
ni; éléments dans le premier ensemble, et de nr; éléments dans le deu- 
xième. Dans les objets anthropologiques, les crânes anciens, par exem- 
ple, il y a souvent des cassures et des déformités, ce qui fait qu'on ne 
peut toujours mesurer toutes les grandeurs de tous ces objets d'un ensem- 
ble; c'est pourquoi les nombres n;; (EIRE ,P) ainsi que les n,;, peuvent 
ne pas être tous égaux. Si x; et y; sont les valeurs moyennes du caractère 
i dans les deux ensembles et si s; est son écart type quand les deux en- 


sembles sont groupés, alors ledit coefficient (CRL = ‘'coefficient of racial 
likeness!') s'écrit 3 
PAR nr .#n x: = y: 
1 121 Û i | 
CGRL == ; 
p 2 ny; Ÿ n); S; 
p 1 2 
“x5 nn D É y; | 
ND, ner S; 
siny, = ny etnrj = n, (i=l,...,p). Les deux faiblesses du CRL sont que, 


premièrement, il ne tient pas compte des corrélations entre les caractè- 
res, et que, deuxièmement, son espérance mathématique dépend fortement 
des nombres d'observations. 


La question de choix entre deux hypothèses se trouve déjà traitée par 
Student quand il a défini son test de signification pour la moyenne d'une loi 
laplacienne. Le problème de tests, avec la considération explicite des hy- 
pothèses alternatives, a connu un grand développement dans les travaux de 
Neyman et E.S. Pearson Li, 18] mais l'insistance a toujours été mise sur 
ce qu'on appelle ''l'hypothèse nulle". En 1936, Mahalanobis [1 1] a proposé 
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une mesure de distance entre deux lois laplaciennes à plusieurs variables 
et avec la même matrice de dispersion. Cette distance, qu'on désigne par 
D?, est une sorte de généralisation de la fonction ‘'t'' de Student, et repré - 
sente en même temps une amélioration du CRL de Karl Pearson. Soient 
((xi: )) et ((yix)) deux ensembles de mesures comme plus haut, Xi, y;i, les 
valeurs moyennes du caractère i dans les deux ensembles respectivement, 
et soit ((a;; )) la matrice pxp de dispersion quand les deux ensembles 
sont réunis. Soit ((a') )}) l'inverse de la matrice ((ai; )). Alors 
D'- 23 ai (%- 5) (7 - 7) 

Si Mi, Vi , Ci (i,j=1,...,p) représentent les moyennes et les éléments 
de la matrice de dispersion (supposée commune) des deux univers lapla- 
ciens à p variables dont les échantillons ont été tirés, et si l'on pose 


2 ij 
À- œ') (ue vi)(u;-v) 
comme la distance entre ces deux univers, alors on montre que l'espéran- 
ce mathématique de D? est, à un facteur près, très peu différente de À 
Ce facteur qui est une fonction de nombres d'observations, est donc in- 
corporé dans la définition de D?. 


L'étude systématique de ce qu'on appelle ‘l'analyse discriminante!'' a 
été commencée par Fisher 5,6,7]en 1936. Il s'est posé le problème sui- 
vant : Soient à notre disposition les deux ensembles ((x;: )) et ((y;x )}) d'ob- 
servations, que l'on sait d'être deux échantillons de deux univers diffé- 
rents; si (x; > ce xp) sont les mesures des p caractères sur un nouvel 
élément, auquel des deux univers cet élément appartient-il? 


Désignons par z; le caractère i. Alors Fisher prend une fonction li- 
2 
néaire 


à coefficients À; arbitraires, et il choisit ces À: de telle manière que z 
soit une fonction qui différencie ie mieux les deux univers, Les z; sont les 
x; dans le premier ensemble, et les y; dans le second. 


Soit W la variation interne de z dans les deux ensembles réunis, et B 
la variation entre les deux moyennes de z. Le rapport B/ W représente 
d'une certaine manière la capacité de z de différencier les deux univers, 
On choisit les À; pour rendre B/W maximum, Dans ses études sur la 
discrimination. Fisher a montré que les travaux antérieurs (comme celui 
de Mahalanobis, par exemple) sont des cas particuliers de cette analyse 
discriminante., 


S'inspirant des méthodes de Neyman et E.S. Pearson, Welch [2 1] a 
traité le problème de discrimination dans sa généralité de la manière sui- 
vante : Etant donné n observations x, ,...,x,, et deux densités de probabi- 
lité f,(x1,...,x,) et f2(x1,...,xn), le problème est de décomposer l'espa- 
ce euclidien à n dimensions en deux régions disjointes, l'une comme ré- 
gion d'acceptation de f,, et l'autre celle de f,. Cette décomposition doit 
être faite suivant un critère donné d'optimalité. Welch a obtenu des métho- 
des de décomposition suivant deux critères différents : l'erreur commune 
minimum, et l'erreur totale minimum. 


Rao [19 ]a généralisé cette analyse au cas de discrimination entre plus 
de deux densités de probabilité. Dans les méthodes de Welch comme dans 
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celles de Rao, on aboutit aux régions d'acceptation définies par des inéga- 
lités sur les rapports de vraisemblance. Or ces rapports constituent des 
résumés exhaustifs pour le choix entre diverses hypothèses, comme l'a 
montré Mourier [16]. 


La définition d'une distance entre deux lois de probabilité n'a pas tou- 
jours été donnée en liaison avec le problème de discrimination proprement 
dit. Le D? de Mahalanobis n'a trouvé sa justification au point de vue dis- 
crimination que dans les travaux postérieurs de Fisher. Depuis quelques 
années, avec Mourier, Rao et Matusita (2, 13] pour les variables aléatoi- 
res réelles, et avec Kullback et Leibler [9] et Komogoroff pour les espa- 
ces abstraits, on voit l'accent mis sur cette liaison entre les distances et 
les erreurs de discrimination, 


L'utilité d'une ‘'distance!! sur les ensembles de lois de probabilité est 
évidente dans les études de convergence de ces lois. Une première défini- 
tion en vue de telles études est celle de Lévy [10] qui prend comme distan- 
ce entre deux lois de répartition de probabilité la plus grande valeur de 
l'intervalle entre les deux courbes dans la direction de la droite x +y =a. 
La question de convergence sera hors de notre présente étude. 


L'introduction des fonctions aléatoires et d'autres éléments aléatoires 
plus abstraits dans la théorie de la probabilité rend nécessaire la généra- 
lisation des méthodes de discrimination et des définitions de distance aux 
espaces abstraits sur lesquels sont définies des mesures de probabilité. 


Dans la présente étude nous nous proposons donc d'examiner la plu- 
part des distances et des ''divergences'' connues jusqu'à présent entre deux 
mesures de probabilité, pour essayer de voir leurs propriétés importan- 
tes et leur interrelations, s'il en existe. Nous nous plaçons dans le cas 
général d'un espace abstrait probabilisé (X,.ÿ# , M ) où M est une classe 
de mesures de probabilité définies sur un corps borélien® de sous- 
ensembles de X. La distance de Lévy ne se prêtant pas à une généralisa- 
tion à un espace abstrait n'a pas été considérée. 


B. PROPRIETES SOUHAITABLES D'UNE DISTANCE 


Une distance définie sur l'ensemble M de mesures de probabilité 
doit, pour être utile dans les applications, être liée à la facilité avec la- 
quelle on peut distinguer deux lois l'une de l'autre, c'est-à-dire que la 
distance doit être en relation avec l'erreur de discrimination tel que plus 
la distance est grande, moins soit l'erreur (voir Mourier [16] }. Comme 
l'erreur de discrimination dépend de la méthode de discrimination, alors 
le choix de la distance en dépendra aussi. 


Il est souhaitable qu'une fonction, pour être une ''distance'', satisfasse 
aux propriétés d'une métrique, et qu'elle prenne toujours des valeurs fi- 
nies. Mais d'après ce que nous venons de dire, la propriété d'être une 
métrique, à elle seule, n'est pas suffisante, Il faut que la distance soit 
une fonction décroissante de l'erreur, et qu'elle prenne sa valeur maxi- 
mum pour deux mesures singulières (Hlv ). 


Une transformation mesurable n'est qu'une extension du concept de 
fonction d'observations. Comme de telles fonctions ne contiennent pas plus 
d'informations que les observations elles-mêmes, elles ne peuvent pas 
diminuer l'erreur. Donc nous postulons qu'une distance sur l'espace des 
mesures ne doit pas augmenter sous les transformations mesurables à 
moins que cette transformation ne soit un résumé exhaustif. 
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Si l'on augmente le nombre des observations on doit être capable de 
distinguer entre deux mesures avec plus de précision. Ceci nous amène à 
un autre critère, à savoir : la distance doit augmenter dans les espaces 
produits de (X,2 ) avec lui-même. 


PRÉLIMINAIRES 

1. Soit X un ensemble d'éléments x de nature quelconque. Nous appel- 
lerons x un point de l'espace X. Si® est une famille de sous-ensembles 
de X telle que 


(siEey , Fey, alors EU Feÿ et En Fev ; 


(ii) si Eeÿ , alors EecŸÿ , où E est l'ensemble complémentaire de E 
par rapport à X, 


alors © est un corps de sous-ensembles de X, Il est évident que si® est 
un corps et si E; (i=1,... ,n) sont n ensembles appartenant à® , alors 


n 
We E; € © 
\= 


Si est telle que l'union de toute suite dénombrable { E;| d'ensembles de 
appartient à , alors est un corps borélien (ou s-corps). Un espace 
X sur lequel est défini un corps borélien ® de sous-ensembies de X sera 
appelé un espace mesurable, et nous le désignerons par le symbole (Rate LV 
On appelle ensemble mesurable tout sous-ensemble de X appartenant à Ÿ . 


2. Une fonction d'ensemble lu définie sur # est une mesure si elle est 
non-négative et complètement additive; est une mesure finie si H(X) < ; 
est une mesure o-finie si l'on peut trouver une suite dénombrable ou fi- 
nie }E;! d'ensembles mesurables telle que Ej = X, et que p(E;)< c 
pour tout i. La mesure u est une mesure de probabilité si u(X) = 1. 


3, Soit (X,® )et(Y, TG) deux espaces mesurables. On désigne par 
X @ Y l'espace produit de X et de Y constitué par tout couple (x,y) où 
xeXetyEe Y. De même pour tout EeŸ et Fev E ©Q F désigne l'en- 
semble de points (x,y) tels que xe Eetye F. Le plus petit corps boré- 
lien sur X @ Y qui contient tous les rectangles est désigné par ® @ T . 
Siu etv sont des mesures définies sur et © respectivement, on obtient 
une mesure m@vsur les ensembles de # Q TG 


4, Soit u,v deux mesures telles que pour tout ensemble E pour lequel 
H(E) = 0, on ait v(E) = 0. Alors v est dite absolument continue par rapport 
àu, et on écrit 

VU 


Si l'on a à la fois u et h <&v, alors on dit que m et V sont équivalen- 
tes, et on écrit H= v. 
Si pour pet V, il existe deux ensemblesdisjoints À etB tels que AUB=X, 


que pour tout ensemble mesurable E les produits (A M E)et(B N E) sont 
mesurables, et que 


H(ANE)=V(BNE):=0, 


alors on dit que Het V sont singulières l'une par rapport à l'autre, Nous 
écrivons vLu. 


5, Une fonction réelle f(x), de point x, est une fonction mesurable (Ÿ ) 
si pour tout nombre réel c, l'ensemble {x nf) e } appartient à ® . 


On sait que 


(Théorème de Radon-Nikodym) si v <= H, alors il existe une fonctionme- 
surable (®Ÿ), f(x), telle que 0 < f(x)<co, et que pour tout E € d, 
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V(E) 6 f(x) du(x) 


Remarques : (i) La fonction f(x) est unique dans ce sens que s'il existe une 
autre fonction g(x), mesurable (Ÿ ), et telle que 


V (E) A g(x) du(x), Eee 


alors 


uix : f(x) + g(x) } En 
on écrit du(x) = f(x) du(x) et aussi f(x) = dv (x)/ du(x). 


(ii) Toute relation suivie du symbole | signifie que cette relation est 
vraie à tout point sauf aux points d'un ensemble de mesure-fi nulle. 


6. Soit MT une classe de mesures de probabilité définies sur à. 
MT est une classe de mesures dominée s'il existe une mesure À (o-finie) 
sur J telle que << À pour toute HE M. La classe M est dite équiva- 
lente à une mesure À ( I = À ) si toutes les mesures dans M sont absolu- 
ment continues par rapport à À} et que À (E) = 0 pour tout ensemble E € 
pour lequel ui (E) = 0 pour chaque Le MC simultanément. 


Toute classe MT dénombrable ou finie est dominée. Halmos et Savage 
[8] ont démontré que si MT est dominée, il existe une mesure de probabi- 
lité À telle que M=X 


Dans ce qui suit nous supposerons que la famille M est toujours do- 
minée et donc qu'il y aura toujours une mesure de probabilité (désignée 
toujours par À) telle que M=X. 


7. Soit (X,®Ÿ )et (Y,G ) deux espaces mesurables dont les points sont 
désignés par x et y respectivement. 


Une transformation T(x) = y de (X,% ) sur (Y, G ) détermineune cor- 
respondance entre les points de X et ceux de Y de telle manière qu'à cha- 
que point de X correspond un seul point de Y, et qu'à chaque point de Y 
correspond au moins un point de X, 


Etant donné un ensemble F£g G, l'ensemble de tous les points de X 
tels que T{(x)€ F est appelé l'image inverse de F, désigné par T''(F). 
La transformation T est une transformation mesurable si pour tout Fe GT 


T'ANLR) ect 
L'ensemble de toutes les images inverses T '(F), où Fe TG, constitue lui 
aussi un corps borélien T''(G) de sous-ensembles de X. On a évidem- 


ment, T'(T)SS 

8, Pour toute fonctiong(y)onpeut définir une fonction gT(x) sur X telle 
que gT(x) = g(Tx). Si g(y) est mesurable (T) alors gT(x) est mesurable 
(T-'(G})). Pour toute mesure u surŸ , on peut définir une mesure UT 
sur © par la relation 


TA (= (Ta RÉSERT 
Les lemmes suivantes sont dues à Halmos et Savage. 


Lemme 1. Sif(x) est une fonction réelle sur X, alors une condition néces- 
saire et suffisante pour qu'il existe une fonction g(y) mesurable (TG) telle 
que f(x) = gT(x) est que f(x) soit mesurable (rte) 

Lemme 2. Sig(y) est une fonction réelle sur Y, et Uk une mesure sur , 
alors 
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Je etrhap ro' (= retenir 


en ce sens que si l'une des deux intégrales existe, l'autre existe aussi et 
les deux sont égales. 


9.Siuest une mesure de probabilité définie sur (X,#), alors toute 
fonction mesurable(Y), f(x), est appelée une variable aléatoire . Soit v 
l'intégrale indéfinie de f(x) par rapport à la mesure u, supposant que f(x) 
est aussi intégrable. C'est-à-dire, 


V(E)=/ fx)ap(x), Eeÿ 
E 
Alors v est une mesure, et VE He 
Soit maintenant T une transformation mesurable de (X,%) sur (Y,T). 
On a 
I bo DT 
Il existe donc une fonction mesurable (T), g(y), telle que 


vT'(r) = / an) du Ty. Fe 
F 


On appelle g(y) l'espérance mathématique conditionnelle de f(x) par rap- 
port à y, et on la désigne par eu (f/y). Si f(x) = X,.(x) est l'indicatrice de 
l'ensemble E ed , alors g(y) est la probabilité conditionnelle de E par 
rapport à y, désignée par Pu (E/ y). 


10, Une transformation mesurable T est appelée un résumé exhaustif 
[3, 4, 8] pour la classe de mesures MC si la probabilité conditionnelle 
Pp, (E/y) d'un ensemble EE par rapport à y est indépendante de H; c'est 
àaldire que si pour tout Eee il existe une fonction p(E/ y) de y telle que 


PLCE/ y) = p(E/y) [ur ; 


et cela pour toute mesure be E M .| Alors nous avons le théorème suivant 
de Halmos et Sauvage [8] : 


Théorème 1. Une condition nécessaire et suffisante pour que la transfor- 
mation T soit unrésumé exhaustif pour une classe YU de mesures (de pro- 
babilité) dominée est qu'il existe une mesure À sur telle que M=X et 
que dw/dA soit mesurable (T=!(T)) pour toute pe NT. 


Soit T un résumé exhaustif, et soit fu) la dérivée de Radon-Nikodym 
d S à 
a. Si l'on désigne par gu(y) l'espérance mathématique conditionnelle de 
fu(x) par rapport à y, alors Halmos et Savage ont démontré que 


69 = ga TO) [x] 
1. “L’'AFFINITE’DE BHATTACHARYA ET LES DISTANCES Y ASSOCIÉES 
Soit MEMC , VEM deux mesures sur (X,® ) et soit 


dp= fu(x)dx, dv=f, (x) dx, 


Pqv= / Ve nee aux) 


Alors la fonction 
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a été introduite par Bhattacharya [1] pour l'étude de la divergence entre 
deux lois de probabilité. Cette fonction mesure, d'une certaine manière, 
la proximité entre deux mesures de probabilité, Il est facile d'énumérer 
les propriétés suivantes de P (uv le 


(i) P(U,v)= P(v,p) 

(ii) 0<Y (p,v}<l 

(iii) P(U,V) = 1 si et seulement si u=V 

(iv) P(H,V) = 0 si et seulement si iv 
En vue de ces propriétés nous appelons P(u,v ), d'après Matusita [12] 
l'l'affinité'' entre pi et v . Une fonction décroissante convenable de P (4, V ) 
peut se suggérer comme une distance entre pi et V. Nous allons en étudier 


quelques-unes. Avant cela, nous allons établir quelques propriétés de 
P(pv ) qui seront utilisées par la suite. 


Lemme 3, Pour tout ensemble E£d , 


PA POECEMCESMONTC 


On a l'inégalité de Holder : 


J'talar.< (ft a (fee 


où f et g sont deux fonctions intégrables quelconques, l1/p + 1/q= 1et 
E ed. 


Mettons f = V £a Ge) 1Q = | f,(x), p = q = 2. On obtient 
REV £ud À)? GAME = Vu(E)v(E) 


Théorème 2, Soit T une transformation mesurable de (X,Ÿ$ ) sur (Y,T), 
et soient gu(y) et g,(y) les fonctions définies par 


duT = gu(y) AXT" 
de Te) d AT 


soit 
P(pv)= P(UTS UT") -[ g(y) (y) x) 
Alors 
PR V) € Pr (pv) 


On peut écrire 


fs COPAEUEUE AT") = /1 SH (7/80) dvT”(y) 
Par la lemme 2, 
F (v)=/ ete Tt) av (x 
) 


où guT(x) et g, T(x) sont mesurables (T-'(T)). Mettons 


69 = ET 


64 B.P. ADHIKARI - D.D. JOSHI 


et soit G{° l'ensemble 


Alors 
Ge (Ch TI 
Soit 
= k 
Sn 237 (ec) 
alors E 


Or, pour xEe GT ;, on a 


Donc 
k? 
2 v (a) €" 0 vb) Er ‘v (ai) 


C'est-à-dire 
k? Lu(Gk) k+1) 


22n v(G) T 21 
ou 

k_ MG) ck+1 

2 V(G)EML A 
On a donc 


x EG auf) = sn = Ne 2 RCE) v (GR) 


Alors, par la lemme 3, 


ne A V£,69 HR) AXE) > J! FOROENS 


u (Go) (Ge) = Pr (uv) 


Halmos etSavage, dans leur démonstration du théorème cité plus haut, ont 
montré que si la transformation T est un résumé exhaustif pour M , alors 
pour toute H 


(x) =e « TR) 
Donc, il est évident que l'affinité P(u,v v) reste invariable sous une trans- 
formation T qui est un résumé exhaustif, 
Une autre propriété de P(u,V) est la suivante : 


Soit (X®@ X,S@J) le produit cartésien de (X,@ ) avec lui-même, et 
A AE V @V les mesures produites correspondant à & et à V |. Alors 
il est facile de voir que 


P(HPV®)= P(HSu, v S)=[P(uN)] 
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et que, de la même manière 


P (at, vin) = ETOI 


Nous allons étudier maintenant les propriétés de trois fonctions de 
P(pv) que l'on peut proposer comme une distance entre Hetv. 


A) d4 (uv) = 1-P(uv) 


1 / [RE - Re] ac 


Cette fonction avait été proposée par Kolmogoroff dans ses conférences à 
l'Institut Henri Poincaré en Novembre 1955. Cette fonction a des proprié- 
tés suivantes : 


1) dy(pv) = d4 (v,H) 
2) o<d{u,v)< i 
3) d,, (p,v) = 0 si et seulement siu=v 
. 4) dy (u,v) =1 si et seulement si ulv 
Celles-ci sont des propriétés d'une ''métrique'' mais l'exemple suivant 
montre que dy(pv) ne satisfait pas à l'inégalité triangulaire. 


Exemple. Soit f, (x), f, (x) deux densités de probabilité laplaciennes 
N(u,,9 , N(b,,02 ), Alors si H= Hz 


PE , 5) = | 29, /(oÿ+ d$ ) 


Prenons trois densités de probabilités laplaciennes 
N(p,c), N(u, 40), Nu, 50) 


Alors 
P (£,, 2) = 0,68 


P(£,, f;) = 0,98 


P (fi, fa) = 0,63 
C'est-à-dire 


PE ED) MP (ES 4) LE PE 63) 
ou 
dy (£,s 1) Fate, Led, ÉTRRUS 
5) D'après le théorème 2, 
1-P(u,v) = 1- P.(H,v) 


C'est-à-dire que dy(lv) n'augmente pas sous les transformations mesu- 
rables et reste invariabie si la transformation T est un résumé exhaustif. 


6) 1 PAU, v®) = 1-[P(pv)] > 1-P(pu,v) 
c'est-à-dire : 
d, (Uvm)=1-[P(piv) > dus v) ; 


si p(p,v)<1, alors 
limd,}(1419 009) 
n —00 
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7) Pour voir la relation entre d,{(u,v) et l'erreur de discrimination, 
considérons la décomposition suivante de l'espace X. 
R =}x E f(x) > (x) 
R =|x : L(x) < f(x) | 


On sait que cette décomposition rend minimum l'erreur totale, c'est-à- 


E£Ÿ 


dire rend minimum l'expression 
u(E) +v(E), 


Or 
P(Hv) = f(x) £,(x) dx, (x) 


3 7 VAE av(x) c 7 AE auto) 


>v(R)+u(R) = € (erreur) 


C'est-à-dire 
Eu 102dà (Ru) 


= 


comme une distance est celle 


B) La deuxième fonction dep{u,v) prise 
proposée par Matusita (2] et s'écrit comme 


d,(uv) aue (EG) -VA Go axp)} 


= RG-p (uv JE 


On vérifie facilement les propriétés suivantes de d, (HV) : 
1) di(u,v) = di(v,u) 


2) o<d(u,v)< V2 
3) di(u,v) = O'si et seulement si p= v. 


4) di(u,v) = V2 si et seulement si uLv 


5) di(u,v) +d;(v,E)>d; (u,.8) 
cette dernière est une conséquence directe de l'inégalité de Minkowski, 
6) di(u,v) n'augmente pas sous les transformations mesurables, et 


reste invariable si la transformation T est un résumé exhaustif. 


D cure) 
feu > af ptites un 


En plus, si p (pv) — |, 
: 3 n) ,ÿ(n) 
lim di(Ht I) )— 1/2 


8) Soit, comme plus haut, € l'erreur totale de discrimination. Alors 


UE j 
7 dj(uwv)=1-p<1-Ee 


11,22 
Œ Æ 1 ---d; (uv) 
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La relation entre di(u,v) et des erreurs dans les tests des hypothèses a 
été étudiée par Matusita et les autres (42% 13514, 15] : 


C) La troisième définition de distance découlant de ‘'l'affinité'' est 


d{u,v) Fri0s p (HV) 


Une forme générale de cette fonction a été étudiée par Chernoff [2] , 
Il propose comme ‘'divergence!!' entre u et v la fonction 


08 | in PACE ax) | 


La distance d,( ,V) satisfait à toutes les propriétés que nous venons 
d'étudier pour les distances précédentes; on doit toutefois faire lesremar- 
ques suivantes : 


1) ddpv) = pour Liv, et reste finie dans les autres cas. 
2) d{H;v) ne satisfait pas l'inégalité triangulaire, 
3) du pl, V0 ) = n dfpiv) 


Il. DISTANCE A L’ERREUR COMMUNE MINIMUM DE RAO 


Soit (E,Ë) une décomposition mesurable de X telle que EUE = X etque 
p (Æ) = v (E) = de 


Si, dans le problème de discrimination entre Lu et V , l'on prend E comme 


région d'acceptation de et E celle de V , alors les deux erreurs de dis- 
crimation sont égales à &.. 


Soit œ= inf œ 


ET : 


Alors Rao Ho] définit la distance entre pu et V comme 
dt) = Water 
On peut énumérer les propriétés suivantes de dr(uv) 


1) d,(uv) = de(v,) 

2) <div) 1 

3) dA(pv) = si et seulement siu=v 

4) dA(pv) = 1] si et seulement si uv 

5) dAlpv) + dA(v,E)>dA(u,E) (voir Rao [19]) 


6° La distance n'augmente pas sous les transformations mesurables. 
En effet, soit 


d EL 
ax ta) Re 


Si (A, A) est une décomposition de l'espace X définie par 


. f(x) RE el 
« > mA ge <k 


AVE 


68 B.P. ADHIKARI - D.D. JOSHI 


où k est choisi tel que 


u (A) = v (A) 
alors l'erreur commune minimum & est donnée par 
= (A) = V (A) 
(voir Welch [21]). 


Soit T une transformation mesurable de (X,®ÿ) sur (Y,T), et soit 
(Ar, Ar) la décomposition de Y définie par 


. 8 (y) 1 AU : gu(y) Le | 
ox | En | rh | 
où | 
du T'= gu(y) dx T, | 


LE 


AvSTI = 8,(y) dv D: | 

et où k!' est choisi tel que | 
= 4, — | 

T'(A;) = vT (AS) | 


Alors l'erreur commune minimum «dans l'espace (Y,T ) est donnée par 
1, “ 
œ=uT (As) =v T(A;) 


Considérons maintenant la décomposition (T' (A,), T '(A,)) de (X,® ). 
On a évidemment 


u(T{A;))=v(T'(A;))= 0e, 


c'est-à-dire que pour cette décomposition les deux erreurs de discrimi- 
nation sont égales, Comme «à est la valeur minimum de cette erreur com- 
mune, alors 

a<X. 


De plus il est évident que la décomposition (T ‘(A,), T '(A,)) de (X,> ) 


rend l'erreur commune minimum si l'ensemble A pour lequel inf x est 
EeY 


atteint est mesurable T_'(T). Dans le cas où T est un résumé exhaustif A 
les régions A,A sont elles aussi mesurables T '(T), puisque les fonctions 
£u(x), f,(x) sont mesurables (T-'(T)). Donc dans ce cas 


X=œ+ 


| 
| 
Ainsi nous voyons que la distance d,(h,v) n'augmente pas sous les | 
transformations mesurables et reste invariable pour les résumés exhaus- | 
tifs. | 


7) Désignons par «!° l'erreur commune minimum de discrimination 
entre les mesures  @u et V@V sur l'espace produit (X® X,#@.7 ). 
Soit R la famille de tous les ensembles Ge Yo pour lesquels 
(2) 


(u ® p) (G)= (V8 v) (G)= a 


Alors, par définition 


2) inf (2) 
at? = GER *G 


Soit À! la famille de tous les ensembles G! de la forme | 
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E ® X 


W(E) =v(E)= 0, 


(HS p) (G')= (VS v)(G)= «= ae 


Alors comme À « R, nous avons 


œa=int œ.= inf CESMinr (West 


œ œ 
Ferme al GER! GER SG 


Donc, si l'on désigne la distance entre u@ u et V@v par CH, VU) 
on à 


due, v®)= 1-01 -x= d{u,v) 


Ainsi on voit que la distance d, augmente avec le nombre d'observations. 


Dans le cas des lois EME à plusieurs variables d AA est 
à la distance D? de Mahalanobis par la relation 


dr(pv) = SN 0 e 2) 


Hi. DISTANCES ASSOCIÉES A ‘ L'INFORMATION” DE SHANNON 


liée à 


Kullback et Leibler [9] ont généralisé l'information de Shannonet ont 
proposé la forme suivante : 


L'information de discrimination de u contre V est 


I (uv) = ee A e x) dx). 
I (V,H) PA) Ex) d'A fx) : 
x 


La ''divergence''" entre H et V est alors définie par Kullback et Leibler 
comme 


De même 


J(u,v) =1(pv) +1(v,p). 


Une définition analogue a été proposée par Mourier Hé] : 


5i-/ [esta] sont [ravi 


Ref fs il £, (69 82,69 [L(v,p)]° 


Alors la distance d{uv) se définit par 


d (pv) _ I(mv) +1 (v 
Zigf > À 


Soit : 


et 
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Nous allons énumérer lespropriétésdeces deux distances successivement. 


(a) La ''divergence!'' J (u,v) de Kullback et Leibler. 

1) 3 (uv) = 3 (0h). 

2)J (pv) >o, et peut atteindre la valeur infinie même pour les 
mesures non- coli ren, 


3) J (p,v) = 0 si et seulement sifi = V. 
4) L'inégalité triangulaire n'est pas toujours satisfaite par J (lv) : 


Exemple : Sif,(x )etf 2 (x) désignent les densités de probabilité . de deux lois 
laplaciennes avec la même moyenne p et les variances g? et ©, alors 


al (of- oi) 2) 


J(1,2)1= Di NT 


Si l'on prend trois lois 
N(u, o), N(u, 40), N(p, 50), 
alors 
JUL,2).=1225/32% J(2,3) = 81/800, (1, 3 57650: 
Donc 
JOUA 2 PAIN 2) (les). 


5)J (,v) n'augmente pas sous les transformations mesurables, mais 
reste invariable si et seulement si T est un résumé exhaustif. 


6) J (u(, V(M) =nJ(u,v). Cette relation n'a évidemment de sens 
que si J (u,v) est finie. 


(b) La distance de Mourier. 

1) d,(u,v) = d,, (v,p). 

2) dy(p,v)>o, et peut être indéterminée. 

3).4}, (pv) = o si et seulement sil=V. 

4) L'inégalité triangulaire n'est pas toujours satisfaite par dv). 
Exemple. Soit 


2x) EN (US Jr, f,(x) AN (Has O2 à f(x) = N (Us O3 } à 
Alors 


Ha) T2 A1) = Te (Ce s?)+ (up) (oi+ c2)] 
Zi ot [te wf,1 (i-0j)"] 
x 2 


9,2 ot 
S 22 0! (Wim), 1 (ot- 02) 
(à + 2 ca O2 


Si l'on prend comme exemple les valeurs suivantes : 


(p= 1, = 1}, (= 10, G;= 3} (H3= 21, (S = 1), 


on voit que 
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d'à.) ve 7 
365 + 12379 
et 
400 
CPE vs 10 
Donc 


dé (152) + 0,0) dy (1,3). 


Nous n'avons pas pu trouver le comportement de dy(u;v) sous les 
transformations mesurables, mais il est facile de voir qu'elle reste inva- 
riable si la transformation T est un résumé exhaustif. La relation entre 
d'y(Hs V ) et l'erreur asymptotique d'un type particulier a été examinée par 
Mourier. 


IV. DISTANCE ‘EN VARIATION” DE KOLMOGOROFF 


Etant donné un ensemble MG de mesures sur (X,©"), on sait que, pour 
deux mesures u etv dans M , la fonction 


dur) = PPPIU(E) = VE) 


est une distance sur M. Si les mesures dans MM sont des mesuresde 
probabilité avec u (X) = 1 pour toute He M , alors on peut écrire 


d, (uv) = Se b(U(E) + V(E)) 


= ER (V(E) - H(E)) 


= NN) = (uv) => |u-v| 


où E',E et [| sont les variations supérieure, inférieure et absolue de 
la mesure &. C'est pourquoi Kolmogoroff appelle d,(u, V ) la distance en 
variation. 


Si est une mesure de probabilité sur (X,Ÿ) telle que HV <= À (on 
peut prendre À comme, par exemple, (pH) Jésalors 
dy(uv) = | RÉMESNTCNSNRE fx) > f(x) | 
2 ÿ Lx à EC) < £,(x) } = H | SE x)  £,(x)| 
LE TATOEOIENCI 
où f(x) = du/ da ; f,(x) = dv/ dx. 
Nous allons énumérer les propriétés suivantes de dy(u:v) : 
tom dede 
2) dy(piv) = o si et seulement siH=vV. 
3) dy(piv) = ] si et seulement si uLv 


4) Soit T une transformation mesurable de (X,2) sur (Y,T); alors 


AGE 
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Donc 
AUUT 0 EEE [HT (F) - VT(F)| 
SUP 
= |HE) -v Œ)] 
arr a IH EE VE] 
= dy (pv) 


C'est-à-dire que d\(p,v ) n'augmente pas sous une transformation mesu- 
rable. 


Soit, maintenant, T un résumé exhaustif. Alors si 


gu(y) = du T'/d\T" 


et 
dd UE /AX TS 


on sait, d'après Halmos et Savage, que 
fu) = gu TG) [A] et f(x) = 8, T(x) [a] 

Donc 

TT -vT"){ y: guy) - e(y)>0} 

u-v)T"!} y : eu(y) - 80} 


= ( 
= ( 
= (u-v){ x: guT(x) - 8y T(x)> o| 
( 


pv)! x : fu x) - f,(x)>0} 

=.d, (pv). 
C'est-à-dire que dy(p,v) reste invariable si la transformation T est un 
résumé exhaustif 


5) Soit UV, deux mesures sur (X, ,Ÿ ), chacune A4» €t Ms Vo 
deux mesures sur (X,,<7% ), chacune € X, . 


Soit du; = f;(x;) dAjet du;= gi(x;) d'A; ,i= 1,2. 
Alors 
1 
Côte MOV) 27 JE Ge) £G)- 86e) 86e) | à (AN) (ur, x) 


X48X2 


XIe x» 
<7/ nee) 6) + 2/]6)-eco| aux) 


F d? (Hs Vs ) Fe K (Les Va À 


Si l'on prend le produit cartésien de (X,Y ) avec lui-même, on obtient 


)-mx)} Abu)+f£6c)-86e)) 86e) é0PAXex) 


DISTANCE - DISCRIMINATION ET RÉSUMÉ EXHAUSTIF 73 


d, (a ; v(® ) < 2 de (lv) 
De même 


durin)end, ty). 


Soit, maintenant, R la classe de tous les ensembles produits de la forme 


E @ X 
où EE , Alors 


Peer 


Donc 


d,(u?,v9)s d,(p@pvev)z Pr (H8)(E)-(v@v) E)| 


>SR |(peu)(E)-(v8v)E)| 
= +9 |BE)-vE)| 
= d, (H.v) 


Ainsi l'on voit que d,(u”, v(") est une fonction positive, non-décrois- 
sante de n, qui est de plus toujours < 1. Donc 


ira ui" D) dl 
n — co 


6) Comme l'avait indiqué M. Kolmogoroff dans ses conférences, la 
distance d,(H,v ) est associée à l'erreur totale minimum de discrimina- 
tion. Soit, en effet, € (1e V ) cette erreur. Alors 


£ (LV) = de (H(E)+v(E)) 
EH {1-((E)-vE)| 
1- SUP (H(E)-v(E)) 


v (Hv) 


7) I1 subsiste des inégalités entre d,(u,v ), de(n,v ) et ''l'affinité!' 
P(uv ). 


(i) Six ete sont l'erreur commune minimum de Welch et l'erreur 
totale minimum respectivement, alors 2 x € . Donc 


de (pv) = 1-x< 1-€/2 => (1 + d,(Hv) ). 
(ii) 


P(H:v) = / fret ren ane) A) RL LAS DATE y @) du (x) 


fa<f 


> ENS me x) <f, (| 
1+ p}x: fu (x) < f, (x) pv x: fut) < fo (x) | 


1-d, (uv) 


(L 


1 
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PROBLÈMES DE DÉPOUILLEMENTS 


par 


P. DUFRESNE 


NOUVELLE EXTENSION DU PROBLÈME DE SCRUTIN 


Probabilité pour que le nombre des bulletins extraits au nom d'un candidat 
soit constamment supérieur à k fois le nombre des bulletins extraits au 
nom d'un autre candidat. 


Nous nous proposons de rechercher la probabilité pour que durant tout 
le dépouillement le nombre des bulletins extraits portant le nom du candidat 
À soit constamment supérieur à k fois le nombre des bulletins extraits 
portant le nom du candidat B, 


Le problème n'a un sens que si k est un nombre positif, nous sous- 
entendrons toujours cette condition, Nous garderons les notations a,b, 
N(a,b) avec les mêmes sens que précédemment, 


Je désignerai par N(a,b)A=>kB le nombre des différents dépouille - 
ments qu'il est possible de constituer avec a bulletins A et b bulletins B et 
dans lesquels le nombre des bulletins extraits au nom de À est toujours 
supérieur à k fois le nombre des bulletins extraits au nom de B 


et par P(a,b)AY> kB la probabilité cherchée, 


Conséquence immédiate des définitions : 


_ N(a,b)AZ>KkB 
P(a,b)AZ KB = Tab) 
ANA 


Lemme : Si a +b=>1 et a kb 
N(a,b)A=kB = Nia-1)A=KkB + N(a,b-l)A>œKB 
La démonstration est identique à celle des lemmes précédents. 


a-kb (a+b)l 


Théorème : Si a>kb et si k est un nombreentier N(a,b)}A=B - Ab “aUbi 


Si a<Ckb ilest clair que N(a,b)A=kB = 0, la formule donnerait un 
résultat négatif donc absurde, 

Si a = kb il estencore clair que N(a,b)A=KkB = 0 la formule proposée 
donne bien un résultat nul. 
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Nous supposons a kb 


La proposition estexacte pour 8 = 1 ; en effet, en vertu de l'hypothèse 
a kb le seul bulletin existant ne peut être qu'un bulletin A, il y a donc un 
dépouillement favorable et un seul, résultat qui serait donné par la formule, 


Nous emploierons un mode de raisonnement par récurrence et nous 
établirons que la proposition est exacte lorsque le nombre des bulletins 
est 0, si elle est exacte lorsque le nombre des bulletins est 8, -1 


En vertu du lemme : 
N(a,b,)A> kB = N(a,- 1,b,)JA=kB + N(a,,b,-1)AZ kB 
et en supposant la proposition exacte lorsque 0 = 6, - 1 


: Max khan (bb el) 
MS ee pes Ml CCR DAT UM 
(du fait que a, kb, il en résulte forcément que a,>kb, + 1) 


a\- kb,+k (a,+b,-1)! 


Neal GR Er En Per TETE 


donc 


kba-tl (abs iabar cekbsr ki (ainb ist 
N PAS RE Zn A 1 1 1 1 1 : 
(2,b,) a bee Re Philo szbre D 1enb/e UN 


ENSECRRPS Al a ; 4 2 = kb, + k b, a,+b.)! 
ab NES ED; api le DES 0 1 


_ (à: - kb a,-a,+(a-kb,)b,+kb, (a; +b.) 
Phi 
\ 


a, + b;-1)(a,+b,) 


_ 2n- kb, (a, + b,) 


a, +betirre, NEA 


ce qu'il fallait démontrer, 


Théorème : Si az kb et si k est un nombre ositif uelc 
a - kb (ab)! ab (asbl PR 
pos TE < N(a,b)A> kB ST PI désignant un nombre entier 


immédiatement inférieur à k. 


Il est clair que N(a,b)JA=KkB<N(a,b)A>jB or, 
nes Jbofatb)l 
N(a,b)=jB = a+b al pb! 
donc 


N(a,b)A> KB? Ryp}e 


Il reste à démontrer que 


a-kb (atb)' SN(a,b)AZ kB 
Ab bi 
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il suffit de le faire pour a >kB car lorsque a = kB, k est obligatoirement 
un nombre entier et le cas a été traité précédemment, 


Nous supposons donc a >kb 
Nous emploierons un raisonnement par récurrence : 


Nous constatons d'abord que la relation est exacte pour la plus petite 
valeur possible de 8 c'est-à-dire pour 8 = 1. En vertu de l'hypothèse a kb 
le seul bulletin existant ne peut être qu'un bulletin À ; il y a donc un 
dépouillement favorable possible etc'estle résultat qui serait donné parla 
formule, 


Nous établirons que la proposition est exacte lorsque le nombre des 
bulletins est 6 , si elle est exacte lorsque ie nombre des bulletins est 
6, - 1. Deux cas sont à considérer : 
1°) a, =kb + 1 (ou a, - 1>kb:) 


en supposant la relation exacte pour 8 = 0, - 1, on obtient : 


a,-kb,+l (a;+b.-1 
N(a,,b-1)A> BD a +b, Ho 


: a,-kb,-1 (a;t+b:-1)l 
OO ES Eng es) nn Cv à 


et comme en vertu du lemme : 


N(a,b,)A= kB = N(a,,b,-1)A=> kB + N(a,-1,b,)A> kB 


| 
N(a,b DAZKkBD ARE (es +b,)! 


ail b! 
suivant un calcul déjà effectué 
2) a,<kb, +1 /(ow a, -11< kb,) 


Il est clair que N(a, - 1,b,)A> kB = 0 


donc en vertu du lemme, 
N(a,,b,)A>kB = N(a,,b,-1)JA=KkB 


et, en supposant la proposition exacte pour 8 = (ue al 


a,-kb,+k (a;+b, -1 
N(a ,b JA kB = N(a b-DAZKB> ET DA 


mais, il a été établi précedemment que : 


2,-kb,+k ere ai-kbi=1 (astb,-1)! _a-kb, (aitb)l 
a, (b,-1)! 


Ta, +b, -1 -1 Lab, -12 (8/2 1)Pb0a,+b: .a,1b, | 
et comme a, - kb, - 1 est négatif 


a,-kb,+lt (a,+b,-1) à a,-kb, (ab)! 
Ta Hell "Oo rb 00 
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d'où il résulte que 


(a:+b.)! 


a:!b.! 


a, - kb, 


> 


N(a,,b,)A=KkB 


UE ane 


Si a >kb si k est un nombre entier 


Théorème : 


Si a>kb et si k est un nombre quelconque (positif) 


a-jb 
a+b 


< P(a,b)A>kB< 


a - kb 
a+b 


ALASÈRE ONE 
mm nm mm mn mm nm | , 
2ES mine S 


Ce tableau donne, en fonction de 


Eu 


(a+b)! 
a! b! 


ne] 
A ON] un 
O N 

(] | al L ll 
Le) q Dour 
u + 28 q : Al 
H n il © nr 
me) nl. 
L o) m il œ 
aq = N a 
> M À k «= 
© — 
À Lu] 7 Il 
o ee A il 

« cd ke FA 
0 = nu] a 
nn. & Ne 
É H 
ss 3 

A a 


J SO 
N 
o | © 
(où 
© a 
Q 
4 
# 
cu 
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Ce tableau donne en fonction de 8 et de b les valeurs de N(a,b)A=——B/4. 
On sait que si a kb 


a-kb (a+b)' a-jb (a+b)! 
a+b Ab N(a,b)a> kB écip (etb}t 


et comme ici k = et j=0 


a-b/4 EH N(a, b)A bre (a+b)! 


a+b a Fra a+tb a!b! 


ainsispour 0 =] 'al=2b; b=T5 


9e le B amd id 
44 618 SN SAT prsr 


199,5<N(6,5)A> B/4<252 


en réalité N(6,5)A>= 246 


Ce tableau donne en fonction de 8 et de b les valeurs de N(a,b)A=—3B-4 
c'est-à-dire les nombres des différents dépouillements qu'il est possible 
de constituer avec a bulletins A et b bulletins B et dans lesquels à tout 
moment le nombre des bulletins A extraits est à quatre unités près cons- 
tamment supérieur à trois fois le nombre des bulletins B extraits. 


Il est facile de démontrer que si a>3b-4 et b différent de zéro 
__a+4-3b (a+b+4)! 
N(a,b)AZ3B-4 = STE 
ainsi pour 0 = 11, a = 8, b = 3. 


gen E6les à 
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LA MÉTHODE DES CASES D'ENTRÉE 


Contrairement à la méthode des erreurs qui permet la recherche de 
nouvelles formules, la méthode des cases d'entrée est seulement un moyen 
de vérification de formules déjà établies censées donner en fonction de 
paramètres, Ô et b par exemple, les nombres des dépouillements satisfai- 
sant toujours telle condition déterminée ou tel ensemble de conditions 
déterminées, 


Encore est-il nécessaire que ces formules ne soient pas quelconques 
mais qu'elles se présentent comme des sommes algébriques de termes de 
la forme 


8-1)! 
(b+k)!(P-b-k-1)] 


! , 
(ou encore de la forme : BK EE s ce qui revient au même puisque : 


g! . p-1)! : 6-1)! 
{b+k)! (8-b-k)! " {b+k)! (0-b-k-1)1 * {b+k-1)! (6-b-k)l 


Nous savons que les nombres de dépouillements qui satisfont toujours 
une condition déterminée ou un ensemble de conditions déterminées peuvent 
être calculés de proche en proche à partir d'un nombre générateur : on 
établit un tableau en fonction de deux paramètres 8 et b par exemple, on 
détermine les cases de ce tableau qui correspondent à des couples de 
valeurs de 0 et de b satisfaisant aux conditions imposées (nous appellerons 
surface intérieure la surface de l'ensemble de ces cases, le reste du 
tableau formant la surface extérieure), on détermine le ou les nombres 
a (qui ne peuvent être que les nombres 1 placés dans les cases 

= 1 1b=O0tett0P= IDE 1) enfin on calcule les nombres des dépouillements 
favorables pour les différentes valeurs possibles de 8 et de b à partir des 
résultats de la ligne précédente. 


Il nous est toujours possible de composer un deuxième tableau donnant 
en fonction des paramètres 0 et b les valeurs d'une expression qui serait 
la somme algébrique de termes de la forme 


or ou de la forme : - 8-1)! 
b+k)! (6-b-k-1)! 7 (b+k)'(0-b-k-1)f 


En effet les valeurs de chacun de ces termes (pour les différentes 
valeurs de 0 et de b) sont celles d'un triangle de Pascal dont le nombre 
générateur est le nombre 1 ou le nombre -1 placé dans la case 0 =1, 
b = -k. Il suffirait donc de reproduire un certain nombre de triangles de 
Pascal dont les nombres générateurs (+1 ou -l) sont tous placés sur la 
ligne 8 = 1 et dans des colonnes dépendant des différentes valeurs de k et 
de faire le total des nombres qui tombent dans les mêmes cases, I sera 
encore plus simple, après avoir inscrit sur la ligne 8 = 1 les nombres 
générateurs des différents triangles de calculer les nombres des autres 
lignes en commençant par celles qui correspondent aux plus petites valeurs 
de 0, un nombre d'une case étant la somme des deux nombres situés sur 


la ligne supérieure l'un dans la même colonne et l'autre dans la colonne 
immédiatement à gauche, 


PROBLÈMES DE DÉPOUILLEMENTS 81 


Sur ce deuxième tableau, nous distinguerons encore une surface inté- 
rieure composée exactement des mêmes cases que celles de la surface 
intérieure du premier tableau et une ou plusieurs surfaces extérieures 
composées des autres cases, 


Il s'agit d'examiner si les nombres placés dans les surfaces intérieures 
des deux tableaux sont identiques, 


Cela se produira lorsque les deuxconditions suivantes seront remplies, 


1°) Les nombres générateurs situés dans les surfaces intérieures des deux 
tableaux sont identiques et placés dans les mêmes cases, 


2°) Aucun nombre de la surface extérieure du deuxième tableau n'entre en 
ligne de compte pour le calcul d'un nombre de la surface intérieure du 
même tableau. 


Revenons sur ces deux conditions : 


Première condition : identité des nombres générateurs situés dans les 
surfaces intérieures, 


Les nombres générateurs du premier tableau sont déterminés directe- 
ment ; trois solutions à envisager selon que pour 0 = 1 les deux dépouille- 
ments possibles (celui qui correspond à b = 0 et celui qui correspond à 
b=1) satisfont aux conditions imposées ou qu'un seul deces dépouillements 
y satisfait. 


Or, le nombre 1 situé dans la case 8 = 1, b = O0 est le nombre généra- 
teur du triangle de Pascal représentatif de l'expression 


6-1)! 
b!(9-b-1)! 
Le nombre 1 situé dans la case 0 = 1, b = 1 est le nombre générateur 
du triangle de Pascal représentatif de l'expression : 


(8-1)! 
(b-1)!(8-b)! 


Conclusion : la première condition est vérifiée : 


a) si les deux dépouillements (0 = 1, b = 0 et 8 = 1, b = 1) sont compatibles 
avec les conditions imposées et que la formule proposée contient à la fois 


6-1)! 6-1 u ; 24 
les deux termes ÊT B-b-1)! et b-1)! -0-b)t ou ce qui revient au même, le 
g! 


terme H1(6-b)T 8-b)] 


b) Si seul le dépouillement 8 = 1, b = 0 est compatible et que la formule 


é ient let (91: 
proposée contient le terme *b!(6-b-1)] 


c) Si seul le dépouillement 8= 1, b = 1 est compatible et que la formule 


=)! 
proposée contient le terme : —S nos 


Deuxième condition : les nombres de la surface extérieure sont sans 
influence sur ceux de la surface intérieure (deuxième tableau). 


Il est facile de déterminer les cases de la surface extérieure dont les 
nombres entrent directement dans le calcul de ceux de la surface intérieure, 
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Ce sont celles qui sont situées sur une ligne supérieure immédiatement 
au-dessus ou immédiatement à gauche des cases appartenant à la surface 
intérieure, 


Nous les appellerons des cases d'entrée - d'où le nom de la méthode - 
car ce sont elles qui donnent accès à la surface intérieure. 


Pour que la deuxième condition soit remplie il importe que les nom- 
bres à inscrire dans ces cases soient systématiquement nuls. 


Nouvelle démonstration de la formule : 


N __(atb)! a+b)! 
ST LE farm)! (b-m)l 


Cette formule a été établie dans les hypothèses suivantes que nous 
conserverons : b>m>1l et a=b-m. 


Nous n'examinerons pas le cas b=m = 1 nous supposerons donc m—1, 


Etablissement du premier tableau, Pour 8 = 1, les deux dépouillements 
possibles : un objet de la famille À ou un objet de la famille B sont favora- 
bles, Il y aura donc deux nombres générateurs : les nombres 1 placés 
danses Cases 0-1 0b =10 (et 1401-21 /2b = 


Aucune difficulté pour tracer la ligne de démarcation entre la surface 
intérieure et la surface extérieure : d'un côté de cette ligne seront les 
cases pour lesquelles a> b-m ou8> 2b-m, de l'autre côté les autres 
cases. 


TABLEAU I 


donnant en fonction de @ et de b 
les valeurs de N(a,b)A > B-4 
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TABLEAU 2 
donnant en fonction de et de b les valeurs 


8! 8! 
de la sommet Bi RAA BUT 


œ 
J 


! 
ms 


1 
CN 
J 
re 


1 
D 
[æ] 


1 
FN 
Oo 


1 
AE 
a 


! 
(+ 
œ 
un 


Etablissement du deuxième tableau : nous reproduisons la ligne de démar- 
cation entre la surface intérieure et la surface extérieure, Nous inscrivons 


! 
les nombres générateurs de {2+b)! 


Ce sont les nombres 1 qui tombentdans lescases 0 = 1, b=0et 8= 1, 


a+b)'! 
1 et de - NB 


(7 
Il 


Ce sont les nombres -1 qui tombent dans les cases 8 = 1, b=m et 
DH b=imtle 


Les nombres générateurs tombant dans les surfaces intérieures des 
deux tableaux sont identiques, la première condition est donc remplie, 


Soient 0, et b, les indices d'une case d'entrée, 
Deux cas sont à considérer : 


- le nombre influencé de la surface intérieure se trouve immédiatement 
en-dessous et ses indices sont 6, + l et b,. 


- le nombre influencé de la surface intérieure se trouve sur la ligne sui- 
vante et dans la colonne de droite et ses indices sont 0\+ 1 et b, + 1. 


mais, ce dernier cas ne peut se produire : en effet 8 et b étant les indices 
d'une case de la surface extérieure 


8,<2b,-m 


et cette inégalité est incompatible avec l'inégalité 


b,+1>>2b,-m +2 
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qui résulterait du fait que 6, + 1 et b +1 seraient les indices d'une case 
de la surface intérieure. 


Reste le premier cas : 
Nous avons toujours 6,<2b, -m 
et en outre 8, +1=>2b, -m 
puisque 8, + 1 et b, sont les indices d'une case de la surface intérieure. 


Pour que ces deux inégalités soient vérifiées il est nécessaire que 
0, =2b,- m ou ce qui revient au même a, = b, - m 


PE ! +b,)! 
Mais sia = b-m la différence ete Rte est nulle 


Les cases d'entrée ne contiennent donc que des nombres nuls ce qui 
vérifie la deuxième condition, 


Nouvelle démonstration de la formule 


AZ B-m 


N(a,b) |B>A-r=-4(@#b)l s{a,b) 4 t(a,b) - u(a,b) + v(a,b) 


Cette formule a été établie dans les hypothèses suivantes que nous 
conserverons : a=b-m, b=a-r, m>1 et r =l1., 


Etablissement du premier tableau, Pour 0 = 1 les deux dépouillements 
possibles : un objet de la famille A ou un objet de la famille B sont favora- 
bles et il y aura deux nombres générateurs : les nombres 1 placés dans 
les cases 0 = 1, b=0 et 8 = 1, b= 1. Il n'y à d'exception que si m= 1 
auquel cas le nombre générateur est le nombre 1 placé dans la case 8 = 1, 
b=0 et si r = l auquel cas le nombre générateur est le nombre 1 placé 
dans la case 8 = 1, b = 1. Si m = r = 1 il n'y a pas de dépouillement favo- 
rable donc aucun nombre générateur, 


Aucune difficulté pour tracer les lignes de démarcation entre la surface 
intérieure et les deux surfaces extérieures ; une ligne marquera la sépa- 
ration entre les cases pour lesquelles a > b-m ou 8 = 2b-m et les autres 
cases ; une autre ligne marquera la séparation entre les cases pour 
lesquelles b>a-r ou Ô=2a-r et les autres cases, La surface intérieure 
sera comprise entre ces deux lignes, les surfaces extérieures de part et 
d'autre, 


Etablissement du deuxième tableau, Nous reproduirons les lignes de 
démarcation entre la surface intérieure et les surfaces extérieures. 
Nous inscrivons les nombres générateurs deLe4b)i ce sont les nombres 
a!b! 
1 qui tombent dans les cases 0 = 1, b=0 et 0 = 1, b=1 ainsi que les 
nombres générateurs de tous les autres termes qui sont de la forme : plus 


ü Dos ere 
o a+h)! (b-h)! 


Ce sont les nombres 1 ou -1 qui sont à inscrire dans les cases : 0= 1, 
b=h et 0=1, b=h+1,. 


! 

Nous remarquerons que les nombres générateurs de er sont tous 
a!b! 

les deux dans la surface intérieure sauf si m = 1 ou le seul nombre géné- 
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rateur tombant dans la surface intérieure est celui de la case 8 = 1, b=0 
et si r = 1 ou le seul nombre générateur tombant dans la surface intérieure 
est celui de la case 0 = 1, b=1, Sin=r=1iln'y a pas de surface 
intérieure, 


Il y a donc identité entre les nombres générateurs tombant dans les 
surfaces intérieures des deux tableaux, la première condition est donc 
remplie, 


s 


Un raisonnement identique à celui qui a été fait pour la démonstration 
précédente montrerait que les cases d'entrée sont celles pour lesquelles 
les indices 0, b, a représentent entre eux une des relations suivantes : 


2b-m 
2a-r 


soit a = b-m qui peut s'écrire 8 
soit b=a-r qui peut s'écrire 8 


Hoi 


! 
Mais si a = b-m la somme eee S(a,b) + T(a,b) - U(a,b) + V(a,b) 
est nulle car : NES 


T(b-m,b) = U(b-m,b) 


= \ 
et S(b-m,b)= RHER + V(b-m,b) 


cette somme, pour des raisons analogues est encore nulle si b = a-r. 


Les cases d'entrée ne contiennent donc que des nombres nuls ce qui 
vérifie la deuxième condition, 


TABLEAU ll 
donnant en fonction 
de 8 et de b 

les valeurs de 
AZ B-4 
B>A-2 


N(a,b) 


CN 


ss = 


[ES 
a 
ot 
ET 
w 
Ne] 
eS 
a 
TS 
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\ Z = 1 39 y = w anod 


(a'e)A + (q'e)n - (a‘e)L + (ae)s - Ten 
ù 
SuuOoSs eT 9p SIN9[EA ST Q 2P 32 ÿ ep UOTJOUOJ ue jueuuop 


z NVATEVIL 
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Nouvelle démonstration de la formule 


E 1 
N(a,b)AZKkB = fab} si a kb et si k est un nombre entier. 


Cette formule peut encore s'écrire si b est différent de zéro, ce que 
nous supposerons : 


_(a+b-1)! el 
N(a,b)A> kB a STD EN 


Etablissement du premier tableau. Pour 0 = 1 un seul dépouillement favo- 
rable possible : il faut que ce soit un objet de la famille A qui soit sorti. 
Donc un seul nombre générateur : le nombre 1 situé dans la case 0 = 1, 
b=0, 


Aucune difficulté pour tracer la ligne brisée de démarcation entre la 
surface intérieure et la surface extérieure : d'un côté de cette ligne seront 
les cases pour lesquelles a=>kb ou0 —(k+1})b, de l'autre côté les autres 
cases, 


Etablissement du deuxième tableau, Nous reproduirons la ligne de démar- 
cation entre les deux surfaces. Nous inscrivons le nombre générateur de 


2: 
2 . c'est le nombre 1 placé dans la case 8 = 1, b = 0 et le nombre 


SH 
générateur de -k Heu c'est le nombre -k placé dans la case 0 = 1, 
bu 1 


TABLEAU ll 


donnant en fonction de 8 et de b les valeurs 
de N(a,b)A >3B 
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TABLEAU 2 


donnant en fonction de 8 et de b les valeurs de la somme 


atb-1)! _, (atb-1)! 
bat li 


Ie 


(] Ü [l 
! 
e] 


1 
ut 
1 
nn 
[æ) 
l 
D 
u1 
! 
1 
ES 
! 
es 


-25 -45 -17 
EE 
1 


a] 


! 
ui 
œ 
[ 
nm IrFeuIRS 
ï 
to 


-23 


| 


Il est clair que seul le premier de ces deux nombres générateurs 
tombe dans la surface intérieure, 


en 


Les nombres générateurs tombant dans les surfaces intérieures des 
deux tableaux sont identiques ; la première condition est donc remplie. 


Soient 8, et b, les indices d'une case d'entrée, 
Deux cas sont à considérer : 


- le nombre influencé de la surface intérieure se trouve immédiatement 
en-dessous et ses indices sont 8, + 1 et b, 


- le nombre influencé de la surface intérieure se trouve sur la ligne sui- 
vante et dans la colonne de droite et ses indices sont 6, + let b, +1 


mais, ce dernier cas ne peut se produire : en effet, 0, et b, étant les 
indices d'une case de la surface extérieure 8,<(k+1)b, et cette inégalité 
est incompatible avec l'inégalité 


8,+1=>(k+1)b, + k 


Reste le premier cas, 
L'inégalité 0,<(k+1)b, subsiste 


et puisque 8, + 1 et b\ sont les indices d'une case intérieure 


8,+ 1=(k+1)b, 
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Pour que les deux inégalités soient simultanément vérifiées, il est 
nécessaire que : 


6,= (k+1)b, ou ce qui revient au même : a, = kb, 


Mais si a = kb la différence 


tn a, +ba-1)! 
a,-1)1B, de Det) est nulle 


Les cases d'entrée ne contiennent que des nombres nuls ce qui vérifie 
la deuxième condition. 
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ANALYCITÉ ET PERIODICITE 
DES FONCTIONS CARACTERISTIQUES 


par 


M. GIRAULT 


Le but de cette note est d'indiquer quelques particularités qui peuvent 
se présenter concernant les propriétés d'analycité et de périodicité des 
fonctions caractéristiques (f.c.). 


Rappelons brièvement les résultats connus de cette question : 


Soit Y, (Z) la f. c. associée à la variable aléatoire X de fonction de 


répartition F (x) 
+00 
g, (2) f e* 4F (x) 


Cette intégrale est absolument convergente pour Z réel et définit alors 


une fonction continue de Z ,. L'intégrale peut exister pour Z complexe = 


Ry . e? dF (x) 


Cette intégrale est convergente en général lorsque y est dans un in- 
tervalle contenant l'origine. 


FÉRIY. 
a)Pourati=A0 


Soit : y'< y < y''. Elle définit une fonction réelle analytique de y. 


b) Si y, est une valeur de l'intervalle de convergence précédemment 


défini : 
itx -y,x 
Lens [2 <a te) 


où e ** dF(x) est la différentielle d'une fonction non décroissante à varia- 

tion bornée, Donc Pa est, par rapport à la variable t une fonction ca- 
170 7 

ractéristique à un facteur positif près et par suite : 


| g (t,yo) < Ÿ (0,Y0) 


Ce résultat implique que l'intégrale qui définit Ÿ{, y) 25t convergente 
dans toute la bande 


Vie L'()< y''où I,,,est la partie imaginaire pure de Z. 
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La fonction ainsi définie est analytique dans cette bande (y'y''). 


J'ai démontré (1) (thèse p. 239) que dans l'intervalle y'y'' de l'axe des 
y la fonction Y (y est convexe : 


Y + Ya 
BAPE VENIR V AVEC — = Y2 
Lo + Los) 
On a : DES 


l'égalité n'étant réalisée que pour WY-=1 


L'ensemble des propriétés précédentes permet d'étudier la périodicité 
de Ÿ (Z) dans la bande de convergence. J'ai démontré (1) que la seule pé- 
riode possible de Ÿ (Z) est une période réelle. Pour que 8 soit période, il 
faut et il suffit que la distribution de probabilité soit une distribution de 


treillis ee , c'est-à-dire que X soit presque certainement égal à 


: 2 
l'une des valeurs : O, ee tee LS SHANE 


A ces résultats ajoutons les remarques suivantes : 


1. REMARQUES SUR LE DOMAINE D’ANALYCITE OÙ BANDE 
DE CONVERGENCE 


Bien entendu, y' ou y'' ou les deux à la fois peuvent être infinis et la 
bande de convergence devient un demi plan ou le plan; elle peut aussi se 
réduire au seul axe réel (c'est le cas de la loi de Cauchy). 


En général, la fonction devient infinie en certains points au moins des 
limites; mais il est intéressant de signaler que Ÿ (Z) peut être bornée 
dans tout le domaine de convergence. 


Exemple : 


C'est bien une f, c. car on lit sur ce développement qu'elle est asso- 


: RUE 1 
ciée à la distribution de masses ——> enx=n 
n 


FA 
où k est défini par D L 3 ter 
n= 2 n k 


Pour Z imaginaire pur, égal à iy, on a : 


= 1 e”? \" 
Pan =K2 . ui 


(1) Les fonctions caractéristiques et leurs transformations - Publications ISUP 
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y 


série convergente pour < 1l.Soit pour y > - log 2 le domaine de 


convergence est le demi-plan I(Z) = - log 2. Or dans tout le domaine y 
compris sa frontière ® (Z) est une fonction bornée. 


Il. REMARQUES SUR LA PERIODICITÉ 


Si $ (Z) ne peut pas avoir de période complexe dans son domaine de 
convergence défini plus haut, il peut toutefois arriver que Ÿ (Z) puisse 
être défini dans l'ensemble du plan complexe sauf en certains points sin- 
guliers ou en certaines lignes singulières et que la fonction étudiée dans 
l'ensemble du plan complexe admette une période complexe, 


1 


Un exemple simple est donné par Ÿ (Z) = 2 


C'est une fonction caractéristique. Elle est associée à une distribution 
continue admettant d'ailleurs pour densité la même expression analytique : 


1 : Prat 
ee or pour Z=iyona Yi, RTE 
7, En = 
2 
la bande de convergence est limitée par y' = - ee et y''= ne 


et la fonction admet la période imaginaire pure Zi. 


EXEMPLE DE FONCTION CARACTERISTIQUE ELLIPTIQUE 


On peut même avoir des fonctions caractéristiques doublement pério- 
diques dans l'ensemble de leur domaine de définition. 


Ainsi soit : 
oo Ter 
Belle pete 2) 


s 


où le produit est étendu à toutes les valeurs impaires de n de - co à + 


Remarquons tout d'abord que : 
LS Lt on 
= a —  —————— 2 
est un produit convergent pour k positif différent de 1 et pour u kan 
(n entier positif ou négatif). 


En effet, supposons par exemple 0 <k 1, pour n = 0 le terme gé- 
néral du produit infini a pour logarithme : 
TELE Ka 


CA 2 RICE n impair = 0 
D+k 1-ku (u-1) nimpair 20 


log 
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et pour n<0 n'=-Î-n 
LE lies 


108 — 7 


nt y -1 c : 
DER D-knut 2k (1) ivre 


La fonction P(u) est une fonction loxodromique qui réalise : 
P (kfu)= P(u) 
La fonction ® (Z) admet donc la période imaginaire pure iy = 4i log k 
mais elle admet aussi la période réelle t = 2T 
® (Z) est une fonction caractéristique analytique dans la bande : 
(y' = - log k y'' = + log k) car c'est le produit des fonctions : 
RS A ls 


AAA Dre EEE 14 


pour toutes les valeurs impaires positives et négatives de n. 


PourEns—=10 


NEA - pn x) 
= ——— ps 
TP A nr (+25 ke 
Et pour n=-n' 0 


AU k'1 = pn! jo 
(= (+25 k e 


Ces développements indiquent les distributions de probabilité corres- 
pondantes, 
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SUR LES INTÉGRALES STOCHASTIQUES 
ASSOCIÉES AUX PROCESSUS PONCTUELS 


par 


RAMAKRISHNAN A. et S. K. SRINIVASAN 


RÉSUMÉ : L'interprétation phénoménologique des intégrales des fonctions 
aléatoires, qui a récemment été formulée par un de nous, trouve son 
application aux problèmes de la Physique où interviennent des processus 
ponctuels plus complexes que ceux du type de Poisson simple. 


INTRODUCTION 


On rencontre en Physique des processus aléatoires que l'on peut 
représenter symboliquement par des intégrales associées à un processus 
plus simple qu'on appelle le processus aléatoire de base (p.a.b.), et dont 
on connait les caractères stochastiques et le mode de développement, Bien 
que l'intégration stochastique ait attiré l'attention des mathématiciens, il 
n'y a pas eu jusqu'ici de tentatives systématiques d'interprétation de ces 
intégrales du point de vue physique. Dans de récents travaux (que nous 
appellerons désormais articles A, B, C), un de nous (R) a développé une 
interprétation phénoménologique de ces intégrales à partir du concept de 
ll'évolution!'' du processus stochastique, en fonction du paramètre qui 
représente la ''durée'' du processus. 


Dans l'article À, en vue de donner une idée des possibilités d'applica- 
tion aux problèmes physiques réels, nous avions choisi le processus de 
Poisson comme p.a.b,., ce dernier étant le processus évolutif non-multi- 
plicatif le plus simple, Ici nous allons étudier des intégrales associées à 
d'autres processus ponctuels, Bien que nous n'insistions pas sur des 
problèmes physiques réels qui correspondent à des intégrales, nos résul- 
tats reposent, toutefois, entièrement sur des raisonnements phénoméno- 
logiques, ce qui leur permet d'être applicables aux problèmes physiques, 
Comme la théorie de l'intégration stochastique se trouve dans un état 
naissant, nous avons considéré utile de présenter nos résultats qui consti- 
tuent un premier pas vers l'application de cette théorie aux problèmes 


physiques concrets, 
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FORMULATION GÉNÉRALE DU PROBLÈME 


Considérons la variable aléatoire yn(t) définie par la formule : 


Ym (t) fe. (tm) 4 taf nain) d'in ‘| &, (ti) nft:) dt: (1) 


où les ®;(t) sont des fonctions certaines connues du paramètre à une dimen- 
sion t ; n(t) désigne le nombre de points aléatoires, ou des évènements, 
dans l'intervalle (0,t) sur l'axe de t. Notre but est de déterminer la densité 
de probabilité (d,p.) de y(t), ou tout au moins ses moments, quand on 
connait les caractères statistiques du processus nt). Dans la partiel 
nous supposerons que n(t) est un processus de Furry, le plus simple 
processus ponctuel multiplicatif dont nous disposons, et ensuite, dans la 
partie Il, nous étendrons ces résultats à des processus ponctuels plus 
généraux, 


PREMIÈRE PARTIE 


1. n (t) REPRÉSENTE UN PROCESSUS DE FÜRRY 


La distribution des évènements aléatoires sur l'axe de t est déterminée 
par la loi suivante : la probabilité qu'un évènement aléatoire arrive entre 
t=T et t=T+dr est {1 +n(tr) }À dr où n(t) est une variable aléatoire 
qui représente le nombre d'évènements aléatoires qui se sont produits 
dans l'intervalle (0,7) sur l'axe de t., La probabilité r(n;t) que n(t)=n 
s'écrit 


mnt) et {lere ts (2) 


Nous supposons qu'il existe un évènement aléatoire à t = 0, qui peut 
s'interpréter comme un individu fl'primaire'' générateur d'une population 
d'individus ''secondaires!'!' correspondant aux évènements aléatoires qui se 
Er dans (0,7). Leur nombre est représenté par nt), et leur d.p. 
Par (2). 


Nous pouvons alors considérer À dt comme étant la probabilité qu'un 
individu en ''produise!! un autre dans l'intervalle (t,t +dt). 


Remarquons que correspondant à une réalisation des évènements à 
t::..., tn dans l'intervalle {9,t), la courbe n°(r) de cette réalisation de 
n(T) se caractérise par des sauts discrets de valeur unité aux points 
tie. tns la courbe restant parallèle à l'axe de t entre deux transitions. 
La valeur réalisée de n(t) s'écrit 

n 
nf{r)=Ù H(r-ti), (o<r<t) (3) 
151 
OÙ t1,...3 tn Sont les points sur l'axe de t où se produisent les évènements 
aléatoires, et H(t) est la fonction de Heaviside, La courbe de réalisation 


de <a peut donc s'interpréter par 


* Conforme à la manière habituelle des auteurs, le même symbole x sera utilisé 


pour désigner toutes les d.p., la distinction entre les diverses d,p. étant évidente selon 
le contexte, 
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Fe Seth. => a (T -t:) (4) 


Nous définissons maintenant y,(t) d'une manière plus générale par la 
formule 


Ym(t) -[a (tm) af à (tm) dtm- nf (to) Dar (5) 


Si bo(to) = 1, alors (5) redevient (1). Renversant l'ordre des intégrales 
SU TO > Directe > Tm». On Obtient : 


mlt) = fl Sato) (ee) dto ; (6) 


où F(t,to) s'écrit comme 


F(t;to) = (6) [a,(e) at [Fe (t2) dt... «[lonttr) dtm (7) 


2. SOLUTION PAR LA MÉTHODE D’ÉQUATIONS AUX DÉRIVÉES 
PARTIELLES 


On voit que l'on ne peut pas obtenir la d.p.n(y-;t) de yn(t) sans recou- 
rir à la d.p. simultanée de l'ensemble des yi(t), i = 0,1,..., m. Pour 
surmonter cette difficulté nous nous servons de la technique développée 
par l'un de nous (R., 1956, article C), et remplaçons t par a dans l'inté- 
grale, a étant un paramètre général quelconque indépendant de t. Si 
n(Ymn;t) désigne la d.p. simultanée de y» et de n, on obtient par un 
raisonnement élémentaire : 


2RGnrnit) 2 x (n+1)n(ymonit) An (Yn -F(a,t),n-l;t) (8) 


avec les conditions initiales 
TN (Ym»130) = 80 8(Y») (9) 
où 02 est la fonction de Kronecker, c'est-à-dire que 06 =0 sin#0,0d6=1, 


et où 6(y,) est la fonction delta de Dirac. Remarquons que dans 
cette équation nous avons affaire avec la fonction de fréquence de 


E (te 
à a to F(a,ts) dt, et non pas celle de y comme définie par (6). En dési- 
[o] Le) 


gnant par G(ym,u;t) et g(s,u;t) les fonctions 
© 


Gymouit) = D r(yminit) u° (10) 
n=0 
= dy 
dons =] GUmuit) à Fed (11) 
[0] 


on obtient 
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dg(suit) - au SE(ssuit); ju e-+Ftb dalsiuit) à gs ste" “D (12) 


Utilisant la méthode habituelle des ''caractéristiques'!', nous obtenons 


= 


la solution à (12) sous la forme 


1 
-Àt RENNES 
gtsrut) = EE af . F(az)-A(t-t) ae (13) 
Q 


u U 


On peut donc obtenir directement la transformée de Laplace* , p(sit), 
den(yn;t). Comme 


7 (Ymi t)= 2 r(Yminit) 
= G(ym,1it) (14) 
on obtient p{s ;t) comme suit : 
p(s:t) = g(s,l;t) 
2 


et FE el e”°F(a;r)-A(t-T) ar (15) 


is 
o 


On remplace a par t pour obtenir la transformée de Laplace dela d.p. 
de ÿym défini en (6). Une méthode alternative pour obtenir la même trans- 
formée qui utilise la distribution simultanée de y» et des n{ti), i=1,...,n, 
est exposée dans l'appendice. 


s 


Les moments de yn(t) se calculent à partir de la relation classique 


E {ya(t) = [(-2Y p(sit)], (16) 


=0 
Les trois premiers moments sont 
Fr 
E { ynlt)} = f et F(t,t)dr (17) 
2 t 
E{y(t)}=2[E{y(t)}] +a f eX F’(t,T) dr (18) 


E{ ele) | = 6 [E{ soft) J+ 6XE{ymtt)] [ @* r'(t,r) ar 
af et Ft,t)dr (19) 


On peut obtenir r(y,;t) explicitement dans les cas particuliers. 
(i) an 1 il ; Po(to)=P:(t) =1 
Dans Ce cas F(t,T)=t.-%" (20) 


Fs Comme dans le cas de x, nous employons le symbole p pour désigner toutes les 
transformées de Laplace, la distinction entre les diverses transformées étant mise en 
évidence par le contexte, 
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Lee) 


NET À_* -(A+s)t}* 
plsst)=eN GAY (1 - e-rat) (21) 
k=0 
x (y, ;t) s'évalue par la formule d'inversion 
1 o+io 
n(yit)=s Fa p(s ;t)e‘)ids 
co 
-Àt Abe N < ne rt) 
sie droit enDe nt DC eue (22) 
k=1 r=0 ? 
x (y, ;t) est continue pour tout y, > O0 sauf Y =ute ny: =trn(y;at) a 


une discontinuité de la grandeur Àe-?Àt 
(ii) m' = 0 ; PE (Th=T, 


Dans ce cas Ft,r) = T ; p(s;t) s'écrit comme 
© 


p(sit)= ns Er; (23) 
k= 


Prenant l'inverse de (23), on obtient 


r (pit) 6(5)+e000 Qt PS BE (1 (SH (y-tr) (24) 


Il est intéressant de remarquer que dans ce cas également, nr (y,it) est 
continue pour tout y, >0 sauf à y, =t ;:n(y.;t) est discontinue à y, =t, la 


=Àt 


discontinuité étant de grandeur Àe Nous remarquons de même que 


E E 
( dn(r) (t - t) dr «| dn(r) TdTt ne représentent pas des processus 
Lr 26bs RAGE: 

lléquivalents'', comme on peut s'y attendre d'ailleurs, puisque la distribu- 
tion des points aléatoires n'est pas ''symétrique'! par rapport aux extré- 
mités de l'intervalle (0,t). Selon la définition de R, deux processus x;(t) 
et x:(t) sont dits ''équivalents!! si les d.p. de x,(t) et de x,(t) sont identiques 
en tout point t. Une discussion plus détaillée des processus équivalents se 
trouve dans les articles A, B et C. D'autre part, quand n(t) représente un 
processus de Poisson, les deux intégrales représentent des processus 
lléquivalents'!', 


3. SOLUTION PAR LA MÉTHODE DES D.P. (!) 


Considérons maintenant la variable aléatoire y,(t) définie par (6), au 
point de vue de la théorie des D.P, (1). Remarquons qu'à chaque point 
aléatoire t; nous associons la fonction F(t,ti). Nous pouvons donc appli- 
quer le théorème démontré par l'un de nous (R, 1953) qui relieles moments 
de Ym aux D.P, des points aléatoires (secondaires), Le moment d'ordre r 
de yh s'écrit : 


(1) Densités-produits. 


100 RAMAKRISHNAN A.. - S.K. SRINIVASAN 


E{ y°(t)] = Al lent F(t,t,) Fltit,) ... Ftite) 


Bibi to dE dt idée 
Ro) sl  F'(t,t:) Fest. . Pitt) AD ter) dt, dt...dt;. 
1 fe 


HD Z chu Î ae) F'(t,t,)...Fht,th)fntast2s e + » stn)dtis dêe » dth 


t1 Jt2 


+ ci [r't,2,) £,(t,) dt, (25) 


La somme sur i pour un h particulier désigne la somme sur toutes les 
diverses combinaisons pour cet h donné, Les coefficients Ci se décompo- 
sent en combinaisons si h > 1 eth<r-1, On écrit Ch(i, où la combinaison i 
appartenant à (h,r) est caractérisée par les nombres 1:,...,1lntels que 
SUP =r} comme 


r Le PR En a l'E EE en ne ee 
CHU Lu Die. till. ll (26) 


où pparmih entiers sont différents l'un de l'autre, et k:,..., kp sont les 
nombres d'entiers d'égales valeurs. 


Il nous reste à obtenir les D.P. fi(t:,..., tHAIE L, ...: T'IpouraIe 
processus de Furry. La fonctionf, (EPA Et) dt, dt est la probabilité 
qu'il y ait un seul point aléatoire dans chaque intervalle (ti, ti + dti), 
i=1l,,,., r. Dans cette définition t:,..,,trne sont pas nécessairement 
ordonnés, De plus, fr est symétrique par rapport à toutes les variables, 
Par des raisonnements élémentaires de la théorie de probabilité on obtient 


tt) = res f(x) dt (27) 
faltisto) = 2 fifti) +a f° fat,T) dr (28)* 


L'équation intégrale pour une forme générale de f; en fonction de 
fi-, peut s'écrire facilement si l'on suppose que t:,..., ti sont ordonnés, 
ce qu'on peut faire sans aucune perte de généralité en vue de la symétrie 
de la fonction par rapport aux variables, 


En supposant t,<t,< ,,..<t,., on obtient 


É(ssbrs on vattheaA Tiriltuta. bi) 
tr 
H | £,(E st +. tre #æ) dt (29) 
o 


Dans la fonction sous l'intégrale, t:,..., tr_1 seuls sont ordonnés, 
T pouvant prendre une valeur quelconque entre 0 et tr. Le premier terme 
de l'équation (29) se déduit du fait que ‘'l'individu!! (point aléatoire) entr 
est produit par un quelconque des r ''individus'' quiconstituent le ''primaire!! 
en t = 0, et par les ‘'secondaires!! en t,,,...,t.:1. Le second membre 
résulte de la remarque que tout individu se trouvant ailleurs sur l'axe det 
peut également en produire un autre en tr, 


Prenant la dérivée par rapport à t- des deux membres en (29), on 
obtient 


of t at AOL E 
rt rt) A TR) (30) 


* Ilest à remarquer qu'à t, = 0 il existe un point aléatoire (primaire) qui est géné- 
rateur de la distribution, et nous n'en tenons pas compte dans la définition des DEP}, 
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dont la solution est la suivante : 
Lettre es Sem astnu) (31) 
Comme fr est une fonction symétrique de t,,,..,t-, alors Q s'écrit 
Q(t. Li de #3" +2 = k(r)e"(t1 +to+:--+tr-1) (32) 


La substitution de (31) et de (32) dans l'équation (29) nous donne la 
relation de récurrence 


kimetrk(r-1) ,; k(T)=À, (33) 
qui admet la solution simple : 
kr} er. | (34) 
Donc 
f(ti,te, ” 4) =Xr| eMHttat:ttr) 
= rl filts) fifte) ... fift,) (35) 


Les moments de y, se calculent par la substitution de (35) dans (25). 
Les trois premiers moments obtenus par cette méthode s'accordent avec 


(17), (18) et (19). 


DEUXIÈME PARTIE 


4. n(t) REPRÉSENTE UN PROCESSUS MULTIPLICATIF GÉNÉRAL 


Keprenons la relation (6) où n(t) représente une distribution des points 
aléatoires pour laquelle la probabilité qu'il se trouve un point aléatoire 
entre t et t+dt est AW(n)dt. La d.p.r(n;t) de n s'écrit 


x (n;t)=W(n-1)W(n-2),..W(0) S A7 er VU (36) 
A%= M {v G) vw (x) (37) 
AN 


Pour obtenir la d.p. de yn(t) comme définie par (6) où a remplace t 
dans F(t,to), nous adopterons la méthode des équations aux dérivées par- 
tielles, Si la d.p. simultanée de y, et de nestn(y,,n;t), par des raisonne- 
ments maintenant familiers, on obtient 


_— »nit) _ -XY (n)r(y, snit)+AV (n-l)n(y,-F(a;t),n-lit) (38) 


Si l'on désigne par p(s,n;t) la transformée de Laplace der(y,,,n:t) 
relative à ym, on trouve 


IAUUE -\Y (n)p(s,n;t) +AY GED et p(s,n-l;t) (39) 


c'est-à-dire 


ë 
p(s,n;it) =AVY CE! p(s,n-1;T) e "At q (40) 
(o} 
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Une solution explicite pour p(s,n;it) n'est possible que dans deux cas 
particuliers : 


(ihseom= 0 et b(t)at 
(ii)les V(n) sont en progression arithmétique 


(i) Dans ce cas (40) devient 


-sT-AY(n)(k-T) 


ë 
p(s,n;t) =ÀY CU p(s,n-l;t)e dt (41) 


Désignant, de plus, par p(s,n;q) la transformée de Laplace de pissnit) 
par rapport à t, on obtient 


p{s,nid) AND PIE mate) (42) 


Nous obtenons, par itération : 


NUS (on 1 E(n=-2)e. (0 


PIE Ca tn-ks +AV(Kk) (43) 


comme 


1 


p (8,059) = ay to) ne 


En décomposant le second membre de (43), nous avons 


, an ohne Eee 45 
1» = AU -l1)Y _ Y (o 
k=0 
où Bk s'écrit 
1 


B= 7 [A{vG)-vG)-(G-0s] (46) 


Prenant l'inverse par rapport à q, 


p(s,nit)=\ V(n-1)W(n-2)...Y(o) SE? exp [-xv (k)t - Ben] 
k=0 47 


Jak 1x 1e 

1 i=r en es 5 k (48) 
= jo I (j-k) Dia G Ch MAR EUAUEE 

rx mr 
où 
j=n (49) 
Gr Ë GE AS A CÉTCIR 

HI 
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Une autre inversion par rapport à s nous donne 


RURAL NEEX V(n-1)W(n-2)...W(o) : c”, 
TI (i-k) Ro 
1x 
EST +0) T0) (,, RE à )] H(P2HF À) (50) 
n(y,;t) est donné par 
nt = D nponit) (51) 


n=0 
(ii) LesWY(n) sont en progression arithmétique 


Soit V(n) -W(n-1) = k. Alors (40) devient 


ë 
p(s,n;t) =\WY (n-1)Y(n-2).. 4(o) e-*Flatns)-AV()(t-tns) 
[o)] 


En-1 ti 
Î e tnt) An lRrtr . ji e”*lato)-aw(1)(t:-t0) p(s,o it )dto (52) 
o 


(] 


x 


où 
p(s,0;to)= e”X Y{o)to (53) 
Donc on peut écrire (52) comme 
E 
p(s,n;t) =\W(n-1)Y(n-2).. V(o)e-xv0x e=ftnpthktni dt, 
Le 
En ki 
Î e”*F(atn-2)+Aktn-2 dénne .f'ersrinerau. dto 
40 o 
€ 
FR b AU PROUE PR AE e"*Y{n)t | [ ef lo,r)+Akr ae| (54) 
D o 
La transformée de Laplace p(s;t) den(y,;t) est 
[ee] 
p(s;t)= ) p(s,nit) (55) 


On peut facilement obtenir les moments de y, à partir de (55). 


4% 
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5. n (t) REPRÉSENTE UN PROCESSUS PONCTUEL DONT LA “ MÉ- 
MOIRE ” NE S’ÉTEND QUE JUSQU'’AU DERNIER POINT RÉA- 
LISÉ 


Considérons la distribution suivante des points sur l'axe de t : la 
probabilité conditionnelle qu'il ne se trouve aucun point entre t =t,et 
t = to(ts > t:), étant donné qu'il y a un point aléatoire en t =t,, estæ(t:-t1), 
et ceci indépendamment de la distribution des points pour t<t,. Il a été 
montré par Ramashrishnan et Mathews (1953) que & détermine d'une 
manière unique la distribution des points sur la droite. 


Remarquons que nue dt représente la probabilité pour que le 


premier point aléatoire se produise entre t et t + dt, étant donné qu'il se 
db(t) 
dt 


trouve un point aléatoire en t = 0. Nous allons désigner - parw(t). 


t 
La valeur réalisée par l'aléatoire Î __ F(t,T) dr est 
o 
n 
F(t,ti) 
at 
pour la réalisation d'évènements aléatoires en ti (i=1,...,n). Le pro- 
cessus d'intégrale n'est pas régénératif, puisque si l'on reporte l'origine 
ent, on obtient pour valeur de réalisation 
n 
F(t-t,,t;-t,) 
i=2 
ce qui ne peut être écrit en fonction de F(t,t;) d'une manière simple, 
= 
Néanmoins, si F(t,ti) n'est fonction que de t-t;*, alors les valeurs 
correspondantes sont 


Ÿ rte) et St) 


et donc le processus est régénératif, puisque les deux quantités ne diffèrent 
que par le terme simple F(t-ti). Dans un pareil cas, l'équation intégrale 


pour n(y,:t) s'écrit facilement, Par la méthode des points régénératifs, on 
obtient 


n (y,it) = fu (T)n (y -F(t-T);t-Tt) dr +[1 [le (Tr) dr ] c(Ym) (56) 


La transformée de Laplace par rapport à y, donne 
L -sF(t-t) $ 
p(s;t) = fu(r)e p(s;t-Tt) dr +[i -[o(n ar] (57) 
(o] 


On obtient les moments de y,, en prenant la dérivée des deux termes 
de (57) par rapport à s. 


* Ceci est vrai si yn satisfait à une équation différentielle linéaire à coefficients 


dn(t) 


constants, > étant le terme non-homogène de l'équation (voir Bartlett, 1955), 
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6. UNE ENTRÉE*x (t) TYPIQUE DU SERVO-MÉCANISME ASSOCIÉ 
A n (t) COMME DÉFINI DANS LA SECTION 5 


Etudions enfin le processus x{t) associé à n(t) de la section 5 de la 
manière suivante 


x(t) = ai tit <ti+: 
ONE PCRe rt. (58) 


où les ai sont des variables aléatoires indépendantes avec la même distri- 
bution X(ai), et où les longueurs d'intervalles ti, - ti = li sont, de même, 
des variables aléatoires indépendantes avec la même d.p.w(l;) d'après la 
section 5, 


Un tel processus stochastique a une grande importance dans l'étude 
des manoeuvres d'un avion dirigé par le radar(James, Nichols et Phillips, 
1947). L'étude de la distribution de 


parts 


est extrêmement difficile et reste encore à faire, En vue des sauts aléatoi- 
res dans x(t) le processus devient non-régénératif même dans lecas spécial 
où F{t,T) n'est fonction que de (t-Tt), ce qui entraine l'échec de la technique 
employée dans la section 5. Nous allons nous borner à la d.p. de y(t) où 


t 
y(t) = [ x(T) dr. La méthode des points de régénération nous donne 
[o] 


n(y;t) DRCUE (y - aift-7);t-T)uo(T) dr 


+ [1 -f'otn ar] ô(y) (59) 


Si p(s ;t) est la transformée de Laplace den(y;t), ona 
t 
p(s 582 ua) J p(s;t-Tr) e 4075  &(r) dr 


+[1- fe (r) ar] (60) 


En principe, (60) détermine p(s ;t). Si nous voulons ies moments de y, 
nous pouvons les obtenir en prenant la dérivée par rapport à s des deux 
membres de (60). Si E|y" fa} SMS) dénlene Re les transformées de Laplace 
de E {yr (t)} et deü{t) respectivement, nous avons 


(y ISE (a) TT (61) 
2 ile 28 6) me 2E{a)elt)id: 
E{ y‘(x) }=E{ a} fe] Te ee E{y ()} (62) 
3 = 3 6w(X) $ 3E{alo(x) d À 
E{ y (Y) }=E{ a fou] 1-0(% dx E{y (9) 
a (63) 


(*) Input. 
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Considérons la variable aléatoire yh(t), comme définie par (6) où n(r) 
représente un processus de Furry. On remarque, d'abord, que si la durée 
de l'intervalle s'étend jusqu'à t, nous associons la fonction F(t,ti) au point 
aléatoire à ti. Donc la variable aléatoire est la somme de ces fonctions. 


La probabilité composée pour que le premier évènement aléatoire se 
“ CE > 
produise entre 7, et T7, + dr, le deuxième entre 7, et Tr + d7>,.,., la nième 
entre Tn et Tn + dt, et qu'aucun autre évènement ne se produise jusqu'à t 
s'écrit 


nl À" exp [-A 7, = 2A(T, -T)-...-nA(T,-7,1)-(n+1)A(t- tn) 


d7, dr, .…. dTn. 


Sir (ymn;t)est la d.p. simultanée de y, et de n, alorsn satisfait à 
l'équation 


t t t 
x(y,,n;t)=\ n! [arf ar. | exp [RA rt -2)(r- 7) >... 
o 4 Tn4 


.….. (nt (t-7)] .6(y% D r(t,n) ) dm (A-1) 


Donc la transformée de Laplace p(s,n;t) den (y.,n;t) s'écrit 
t t t n 
p(s,n;t)=}1" nl [arf de ". exp [-A(n+1} +2 AT- sF(t,r;) ler, (A-2) 
1 n 3 


c'est-à-dire 


pis,nsth= 0" [A f'exp | -X(t-T) - sF(t,t)}ar] (A-3) 


La transformée de Laplace p(s ;t) de x (y,;t) est donnée par : 
@ 


D ptsonit) 


en [2 - x [exp -À (t-T) - s(t,r)jae] (A-4) 


p(s;t) 


[l 


ce qui concorde bien avec l'expression (15), 


ANALYSES D'ARTICLES 
PARUS DANS DES REVUES ÉTRANGÈRES 


JOURNAL OF THE ROYAL STATISTICAL SOCIETY 
SERIES B - METHODOLOGICAL 


Volume 18 - N°1 - 1956 


SOME TESTS OF SIGNIFICANCE WITH ORDERED VARIABLES par F.N, DAVID 
and N.L. JOHNSON. 


Les auteurs indiquent quelques tests de signification relatifs à la loi normale 
(moyenne, écart-type, comparaison de moyennes et d'écarts-type ; test de norma- 
lité), applicables à des échantillons ''censurés'! de deux types 


Type À :Ilest décidé d'arrêter l'expérience lorsque dans un échantillon donné, 
le pourcentage des éléments remplissant certaines conditions a atteint une valeur 
fixée. 

Type B : La durée de l'expérience est déterminée à l'avance. Des fonctions 
des variables échantillons ainsi obtenues et ordonnées, permettent de tester les 
hypothèses étudiées. 


A SEQUENTIAL TEST FOR RANDOMNESS OF INTERVALS par D.J. BARTHOLOMEW 


Cet article étend à une procédure progressive, celle étudiée par MORAN en 
1951, relativement à la division aléatoire d'un intervalle. 


En particulier, lorsque des évènements surviennent au hasard dans le temps, 
il s'agit de tester l'hypothèse que les intervalles entre deux évènements consécutifs 
ont une distribution exponentielle contre une distribution de Pearson de type III. La 
méthode et les résultats de Wald servent à déterminer complètement le test, 


ON THE ESTIMATION OF SMALL FREQUENCIES IN CONTENGENCY TABLES par 
FT GOOD' 


Pourestimer les fréquences correspondant aux entrées des Tables de Contin- 
gence, des difficultés surgissent pour la détermination des fréquences voisines de 
0. L'auteur définit un ''facteur association!' entre lignes et colonnes, pour chaque 
cellule ; ce facteur est égal au produit de la fréquence attribuée à la cellule étudiée 
par celles relatives à la ligne et colonne qui déterminent cette cellule. Deux dis- 
tributions de ce facteur sont considérées : une distribution de Pearson de type III 
etlogarithmico-normale. La valeur moyenne de la même puissance du facteur d'as- 
sociation multiplié par la puissance n° du logarithme est calculée. 


108 JOURNAL OF THE ROYAL STATISTICAL SOCIETY 


ALGEBRIC THEORY OF THE COMPUTING ROUTINE FOR TESTS OF SIGNIFICAN- 
CE ON THE DIMENSIONALITY OF NORMAL MULTIVARIATE SYSTEMS par F.E. 
BINET and G.S. WATSON. 


Avant de classer des individus dans plusieurs populations à plusieurs variables 
ilest coutume de faire des tests préliminaires de signification relatifs aux relations 
structurales qui pourraient exister entre ces populations. FISHER et WILKS ont 
donné deux méthodes. Les auteurs reprennent les calculs théoriques concernant ces 
procédés afin de les comparer. 


SOME NOTES ON ORDERED RANDOM INTERVALS par D. BARTON and F. N. 
DAVID. 


Les auteurs reprennentles méthodes de Whitworth pour étudier la probabilité que 
(n - 1) points étant rnarqués sur un segment unité, r segments aient une longueur 
inférieure à une quantité fixée. En particulier Whitworth utilise des coefficients 
œ...æ) dont les auteurs calculent les moments et coefficients de corrélation. Le 
calcul de la probabilité cherchée est effectué à l'aide de la fonction génératrice des 
moments, ce qui permet celui des moments, coefficients de corrélation, détermi- 
nation de la rapidité de convergence asymptotique etc... ainsi que la distribution 
du rapport et de la différence de deux intervalles ordonnés. 


À TEST OF SIGNIFICANCE FOR AN UNIDENTIFIABLE RELATION par P.A. P. 
MORAN. 


Lorsque deux variables v et w sont liées linéairement avec des coefficients 
inconnus, le problème est de tester l'hypothèse que la relation linéaire passe par 
un point donné du plan ; soit w = &« v contre w = œv +7. 


Ce test utilise des valeurs expérimentales sujettes à des erreurs et la posi- 
tion du point considéré par rapport aux deux lignes de régression empiriques. 


CONFIDENCE LIMITS FOR THE GRADIENT IN THE LINEAR FUNCTIONAL RELA- 
TIONSHIP par M.A. CREASY. 


Cet article fait suite à celui de LINDLEY 1947 relatif aux relations structu- 
rales linéaires lorsque les variables observées sont sujettes à des erreurs. Des 
estimations de l'angle qui détermine la direction de la relation structurale sont 
étudiées ainsi que les intervalles de confiance correspondants. 


ON A TEST OF SIGNIFICANCE IN PEARSON'S BIOMETRIKA TABLES n° 11 par 
S.R., FISHER. 


Le problème de BEHRENS relatif au test d'égalité des moyennes de deux po- 
pulations normales dont les variances sont inconnues, fait intervenir le rapport de 
ces variances, Lors du calcul des Tables de Pearson, WELCH a proposé une solu- 
tion dont FISHER montre quelques défauts qui entrainent des inexactitudes dans les 
tables de Pearson. (FISHER envisagerait l'utilisation d'une information àpriori sur 
0e js 
œ 
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ECONOMIC CHOICE OF THE AMOUNT OF EXPERIMENTATION par P.M. GRUNDY 
M.J.R. HEALY and D.M, REES. 


Lorsque l'on veut étudier l'augmentation d'une production résultant d'un pro- 

cédé nouveau, il est nécessaire de déterminer l'étendue des essais préliminaires , 

afin d'en minimiser le coût, Au lieu de réaliser une seule expérience à grande 

échelle, les auteurs envisagent une succession progressive d'expériences, En rai- 

de la difficulté des calculs, le problème est limité à deux étapes : une première 

suite de n,essais est réalisée et le coût est k' n, (k' connu) ; si une décision ne peut 
N 


être prise, le nombre n, optimum des essais supplémentaires à envisager avec un 
coût k est calculé en fonction de k!', k,n. 


STATISTICH VIERTELJAHRISSCHUFT (VIENNE) 


L. SCHEMETTERER 


1.- Règle de Bayes et méthode de vraisemblance 1951 vol. 4, 3/4 
2.- Sur la règle de Bayes 1952%ol.5; 4 


La plupart des auteurs prennent position pour le théorème de Bayes ou pour 
la méthode de la vraisemblance. Pourtant les deux méthodes ne diffèrent pas es- 
sentiellement. Les divergences disparaissent d'ailleurs si l'on se place du pointde 
vue purement mathématique et n'apparaissent que dans les applications et les inter- 
prétations. 


I1 est bien connu que l'on reproche à la règle de Bayes l'emploi des probabi- 
lités à priori essentiellement inconnues. Invoquant le principe de raison suffisante 
on suppose ces probabilités à priori réparties uniformément. Il y a aussi une rai- 
son mathématique, à savoir que l'emploi d'une distribution à priori différente 
entrainerait des conséquences fâcheuses pour la détermination des probabilités à 


posteriori. 

Pourtant l'hypothèse d'une répartition uniforme des probabilités à priori n'est 
pas suffisamment fondée. L'auteur montre qu'il existe une famille de distributions 
à priori comportant comme cas particulier la distribution uniforme et quin'altère 
pas la structure de la fonction donnant la distribution à posteriori. 

En effet, soit une population à une caractéristique A, p la proportion (incon- 
nue) des éléments possédant cette caractéristique A, l-p la proportion des éléments 


possédant la caractéristique contraire À, Au moyen de n tirages successifs on a 
constitué un échantillon de k élements A et de { éléments A.(k + À = n) La densité de 


probabilité à postériori de p est alors donnée par : 
d(p)rk(1-p)° 

Ja (p) p*(1 - p)° ép 

d(p) étant la densité de probabilité à priori de p (0 < p < 1) 


Y(p;k,l)= 


Dans l'hypothèse de l'uniformité, d(p) = 1, et l'on obtient : 
p<(l- p}° 


Y,(p;k, 1) ER URR 
Jr* (1 - pt dp 
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En choisissant : 


k (l 
mp (cp) 
d (PES Ge rasha) 


il vient : 
k+ky (1 = p}êtls 
ik, ?) = Cp - Y(Pik+k,l+b) 
Jr - p)tt4 &p 


c'est-à-dire la structure de la fonction W,(p;k,{)correspondantaucas oùladistribu- 
tion uniforme des probabilités à priori est inaltérée. 


E, MICHALUP - Sur la formule de Stirling concernant les factorielles 1949, vol.2, 
314, 


On connait la formule du développement de n| 


(Ci). ME RE nes Br _, __Ba sPé | 


+ 
l2#n 083 4n 055.06 .n8 


ve B4 = - TL etc. Le premier 


terme de ce développement constitue la formule bien connue de Stirling : 


DEN V2/Tobe ET 


L'auteur propose le développement suivant 


. 5 n 265n 
à as n on A —, 
(@) nl=V/2rne" n"exp É (30 n2+1) 14 (30 Be 


Si l'on ne prend que le premier terme entre crochets, on obtient une valeur 
plus exacte que celle que l'on obtient en prenant les deux premiers termes du dé- 
veloppement (1). 


où les Bj désignent les nombres de Bernouilli : B2 = 


# = 


PRES 


à Se: 


U 


( 
our 

14 | 
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